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(N ' 

j^jy We derive first-order relativistic dissipative hydrodynamic equations from relativistic 

r—{ , Boltzmann equation on the basis of the renormalization-group (RG) method. We introduce 

a macroscopic-frame vector (MFV), which does not necessarily coincide with the flow velocity, 
to specify the local rest frame on which the macroscopic dynamics is described. The five 
[~^. | hydrodynamic modes are naturally identified with the same number of the zero modes of the 

linearized collision operator, i.e., the collision invariants. After defining the inner product in 
the function space spanned by the distribution function, the higher-order terms, which give 
i-Q ' rise to the dissipative effects, are constructed so that they are precisely orthogonal to the 

zero modes in terms of the inner product: Here, any ansatz's, such as the so-called conditions 
of fit used in the standard methods in an ad-hoc way, are not necessary. We elucidate that 
'■ the Burnett term dose not affect the hydrodynamic equations owing to the very nature of 

the hydrodynamic modes as the zero modes. Then, applying the RG equation, we obtain the 
hydrodynamic equation in a generic frame specified by the MFV, as the coarse-grained and 
covariant equation. Our generic hydrodynamic equation reduces to hydrodynamic equations 
in various local rest frames, including the energy and particle frames with a choice of the 
MFV. We find that our equation in the energy frame coincides with that of Landau and 
on : Lifshitz, while the derived equation in the particle frame is slightly different from that of 

Eckart, owing to the presence of the dissipative internal energy. We prove that the Eckart 
equation can not be compatible with the underlying relativistic Boltzmann equation. The 
proof is made on the basis of the observation that the orthogonality condition to the zero 
qq . modes coincides with the ansatz's posed on the dissipative parts of the energy-momentum 

■ tensor and the particle current in the phenomenological equations. We also present an 

analytic proof that all of our equations have a stable equilibrium state owing to the positive 
definiteness of the inner product. 
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§1. Introduction 

Relativistic hydrodynamic equation is widely used in various fields of physics, 
especially in high-energy nuclear physics 1 ) -3 ) and astrophysics, 4 )~ 6 ) and it seems that 
the study of the relativistic hydrodynamic equation with dissipative effects is now 
becoming a central interest in these fields. 

For instance, it was shown that the dynamical evolution of the hot and/or dense 
QCD matter produced in the Relativistic Heavy Ion Collider (RHIC) experiments 
can be well described by the relativistic hydrodynamic simulations. 2 )' 3 ) The sugges- 
tion that the created matter may have only a tiny viscosity prompted an interest in 
the origin of the viscosity in the created matter to be described by the relativistic 
quantum field theory and also the dissipative hydrodynamic equations. Also the 
relativistic dissipative hydrodynamic equation has been applied to the various high- 
energy astrophysical phenomena, e.g., the accelerated expansion of the universe by 
bulk viscosity of dark matter and/or dark energy. 5 )' 6 ) 
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It is, however, noteworthy that we have not necessarily reached a full understand- 
ing of the theory of relativistic hydrodynamics for viscous fluids, although there have 
been many important studies since Eckart's pioneering work. 11 ) 

1.1. Fundamental problems with relativistic hydrodynamic equation for a viscous 
fluid 

We may summarize the fundamental problems on the relativistic dissipative 
hydrodynamic equations as follows: (a) ambiguities and ad-hoc ansatz's in the def- 
inition of the flow velocity; 8 )' 10 )~ 16 ) (b) unphysical instabilities of the equilibrium 
state; 17 ) (c) lack of causality. 16 )' 18 )- 20 ) 

One might consider that the definition of the flow velocity should be merely a 
kind of the choice of the coordinate space to describe dynamics in a easy and practical 
way, and one can go back and forth between two different definitions by a Lorentz 
transformation. Furthermore, one might consider that the unphysical instabilities 
of the equilibrium state should be attributable to the lack of causality, and one 
can restore the instabilities automatically by solving the causality problem. Thus, 
although the relativistic dissipative hydrodynamic equation has been attracting a 
great interest, it seems that most works are concerned with the causality problem 
and examine phenomenological or semi-phenomenological causal equations with some 
foundations. We argue, however, that the first two problems, (a) and (b), and the 
third one, (c), have different origins, and the first two must be resolved before the 
third one is addressed. The present paper is concerned with the first two problems, 
and the third one will be studied in the forthcoming paper. 21 ) 

(a) The form of the relativistic dissipative hydrodynamic equation depends on 
the definition of the flow velocity, which is equivalent to the choice of the local rest 
frame of the fluid. Typical rest frames include the particle frame and the energy 
frame, and a phenomenological equation for the respective frame is constructed by 
Eckart 11 ) and Landau and Lifshitz, 12 ) respectively. Here, the phenomenological con- 
struction is based on the following three ingredients; (A) the particle-number and 
energy-momentum conservation laws, (B) the law of the increase in entropy, and 
(C) some specific assumptions on the choice of the flow. The first and second points 
are reasonable and used also in the construction of the non-relativistic Navier-Stokes 
equation, while the third point is specific for the relativistic case. To be more explicit, 
let 8T^ U and 5N^ be the dissipative part of the energy-momentum tensor and the 
particle current, respectively. The point is that the forms of bT^ v and 5N^ are not de- 
termined uniquely only by the particle-number and energy-momentum conservation 
laws and the law of the increase in entropy without some physical ansatz's involving 
the flow velocity with u^u^ = u^u v = 1 and g^ v = diag(+l, —1, —1, — 1): 
The implicit assumptions made by Eckart 11 ) are 

(i) be = UfM ST^ u v = 0, (1-1) 

(ii) Sn = SN* 1 = 0, (1-2) 
(m) u^ = A^5N u = 0, (1-3) 

where A^ u = g^ v — u^u u . The condition (i) claims the absence of the internal energy 
be of the dissipative origin, while (ii) and (iii) require no dissipative particle-number 
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density 6n and current z/^, respectively. On the other hand, those of Landau and 
Lifshitz 12 ) consist of (i), (ii), and 



In physical terms, (i) and (iv) claim the absence of the internal energy be and current 
of dissipative origin. As one sees, the first two ansatz's are common for the two 
equations, and the third ones ((hi) and (iv)) are supposed to specify the respective 
local rest frame of the flow velocity. We note that the conditions (iii) and (iv) and 
hence the two of frames specified by these conditions can not be connected with each 
other by a Lorentz transformation. It is here noteworthy that there is a proposal by 
Stewart 14 ) for the condition for the particle frame, as given by (ii), (iii), and 



where Sp = —A^ ST^ V /3 is the dissipative pressure to be identified with the standard 
bulk pressure. Here, the condition (i) of Eckart is replaced by the different one 
(v), which claims a constraint between the dissipative internal energy Se and the 
dissipative pressure Sp. One may ask if both the Eckart and Stewart ansatz's make 
sense or not. It is noteworthy that the most general derivation of the hydrodynamic 
equation on the basis of the phenomenological argument gives a class of equations 
which can allow the existence of the dissipative internal energy Se and the dissipative 
particle-number density Sn as well as the standard dissipative pressure dp, as shown 
by the present authors; 22 ) a brief recapitulation of Ref.22) is given in Appendix A. 

(b) The unphysical instabilities of the equilibrium state might be attributed to 
the lack of causality, and Israel-Stewart's formalism is presently being examined in 
connection to this problem. 8 )~ 10 )' 20 ) Although their equation may get rid of the 
instability problem with a choice of the relaxation times, as shown in Ref.17), we 
emphasize that there exists no connection between the unphysical instabilities and 
the lack of causality. In fact, the Landau-Lifshitz equation is free from the instabil- 
ities of the equilibrium state in contrast of the Eckart equation. Furthermore, one 
should notice that the causal equation by Israel and Stewart is an extended version 
of the Eckart equation and hence it can naturally exhibit unphysical instabilities 
depending on the values of transport coefficients and relaxation times contained in 
the equation. 23 ^ 

1.2. Strategy of the present work: separation of the scales of dynamics and con- 
struction of hydrodynamic equation as the infrared effective dynamics 

Here, we note the hierarchy of the dynamics of the time evolution of a many- 
body system: In the beginning of the time evolution of a prepared state, the whole 
dynamical evolution of the system will be governed by Hamiltonian dynamics that 
is time-reversal invariant. When the system becomes old, the dynamics is relaxed 
to the kinetic regime, where the time-evolution system is well described by a trun- 
cation of the BBGKY (Bogoliubov-Born-Green-Kirkwood-Yvon) hierarchy; 24 ) the 
Boltzmann equation composed of one-body distribution function describes a coarse 
grained slower dynamics, in which time-reversal invariance is lost. Then, as the 



(iv) Q^ = A^5T»Pu p = 0. 



(1-4) 



(v) ST» =Se-3Sp = 0, 



(1-5) 
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system is further relaxed, the time evolution will be described in terms of the hydro- 
dynamic quantities, i.e., the flow velocity, the particle-number density, and the local 
temperature. In this sense, the hydrodynamics is the infrared asymptotic dynamics 
of the kinetic equation. 

Now the relativistic Boltzmann equation is a typical kinetic equation, which is 
manifestly Lorentz invariant and free from the instability and causal problems. 18 ) 
Thus, one sees that a natural way to resolve the ambiguities in the definition of the 
flow velocity and the unphysical instabilities of the equilibrium state is to derive 
the relativistic dissipative hydrodynamic equation from an underlying relativistic 
kinetic equation. Indeed, there have been some vigorous attempts to derive the phe- 
nomenological equations from the relativistic Boltzmann equation; for instance, with 
use of the Chapman-Enskog expansion method 25 -* and the Maxwell- Gr ad moment 
method. 26 ) We would say, however, that these works in the microscopic approaches 
are not fully satisfactory: Although the past works certainly succeeded in identify- 
ing the assumptions and/or approximations to reproduce the known hydrodynamic 
equations by Eckart, Landau and Lifshitz, Stewart, and Israel, the physical mean- 
ing and foundation of these assumptions/approximations remain obscure, and thus 
the uniqueness of those hydrodynamic equations has never been elucidated as the 
long-wavelength and low-frequency limit of the underlying dynamics. Indeed, the 
standard derivation of relativistic hydrodynamic equations based on the Chapman- 
Enskog expansion or Maxwell-Grad moment method 18 ) utilizes the ansatz's given by 
(i) ~ (v) as the constraints on the distribution function as the solution of the rela- 
tivistic Boltzmann equation, rather than consequences of the derivation. Their va- 
lidity or the fundamental compatibility with the underlying Boltzmann equation has 
never been questioned nor addressed. This unsatisfactory situation rather reveals the 
incompleteness of the Chapman-Enskog expansion method and the Maxwell-Grad 
moment methods themselves as a reduction theory of the dynamics. 

Thus, the form of relativistic dissipative hydrodynamic equations is still contro- 
versial and far from being established. The origins of the difficulty of the derivation 
are identified as the absence of an appropriate coarse-graining method that keeps 
the Lorentz covariance and is applicable to the relativistic Boltzmann equation. It 
should be mentioned here that van Kampen 16 ' applied his reduction theory to derive 
a relativistic hydrodynamic equation. The resultant equation was, unfortunately, of 
a noncovariant form. 

In this work, we try to derive the relativistic dissipative hydrodynamic equa- 
tion from the relativistic Boltzmann equation in a more natural and systematic way: 
For that, it is essential to adopt a powerful reduction theory of the dynamics. 27 ) 
As such a reduction theory, we take the "renormalization-group (RG) method". 
In Ref.'s 28)-35), it has been shown that the RG method is a powerful resumma- 
tion method and also gives a systematic reduction theory of the dynamics leading 
to the coarse graining of temporal and spatial scales, which are the key concepts 
in the construction of infrared effective theories *). Indeed, the RG method is al- 



*' A brief account of the RG method is given in Appendix B using an example for self- 
containedness. 
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ready applied satisfactorily to the derivation of the Navier-Stokes equation from the 
(non-relativistic) Boltzmann equation, with no heuristic assumption. 34 )' 35 ) The RG 
method thus should be most suitable for the present purpose to derive covariant 
relativistic dissipative hydrodynamic equations. It is also expected that the physical 
meanings and the validity of the ansatz's posed in the phenomenological derivation 
will be elucidated in the process of the reduction in the RG method. 

In Ref.36), the present authors and K. Ohnishi applied the RG method to derive 
the relativistic dissipative hydrodynamic equations for the first time: A macroscopic- 
frame vector was introduced with which the derivation of a coarse-grained covariant 
equation is made possible, and thus the so-called first-order relativistic dissipative 
hydrodynamic equations were successfully derived, in which the dissipative effects are 
taken into account up to the first order. The five hydrodynamic variables naturally 
correspond to the zero modes of the linearized collision operator. The deviations 
from the local equilibrium are given by the functions that are precisely orthogonal 
to the zero modes with an inner product for the distribution functions. It was found 
that the various local rest frames of the flow velocity can be realized by a choice of 
the macroscopic- frame vector. Subsequently, the present authors 37 ) showed that the 
derived hydrodynamic equation in the particle frame is different from the Eckart and 
Stewart equations, and has the stable equilibrium state in a numerical calculation 
for a rarefied gas where dissipative effects are most significant. 

1.3. Purpose of the present paper 

The purpose of this paper is as follows: (1) We present a detailed and full ac- 
count of the derivation of the first-order relativistic dissipative hydrodynamic equa- 
tions in generic frames from the relativistic Boltzmann equation on the basis of 
the RG method. Moreover, (2) we elaborate and revise some parts of the deriva- 
tion and thereby make it more transparent: We clarify the essential importance 
to properly define the inner product of the distribution functions, in particular to 
make it positive-definite, which was not recognized in the previous presentation. We 
also show using simple properties of the inner product that the usually problem- 
atic Burnett term does not contribute to the hydrodynamic equations on account 
of the fact that the hydrodynamic modes are described by the zero modes of the 
linearized collision operator. (3) We show that the ansatz's (i) ~ (v) for the dissi- 
pative currents exactly correspond to the orthogonality conditions of the perturbed 
distribution function to the zeroth-order solution, and that the Eckart constraint (i) 
~ (iii) can not be compatible with the underlying Boltzmann equation as a corol- 
lary. (4) We also fully examine the properties and its advantageous nature of the 
derived equation in a generic frame. We give a general proof without recourse to 
numerical calculations that the hydrodynamic equations obtained in our formalism 
have the stable equilibrium state even in the particle (Eckart) frame on the basis of 
the positive-definiteness of the inner product. 

This paper is organized as follows: In §2, after a brief account of the basic 
ingredients of the relativistic Boltzmann equation, we introduce the macroscopic- 
frame vector, and summarize the ad-hoc aspects in the standard methods such as 
the Chapman-Enskog expansion and Maxwell-Grad moment methods. In §3, with 
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use of the RG method, we reduce the relativistic Boltzmann equation to a generic 
form of the relativistic dissipative hydrodynamic equation, whose frame is not speci- 
fied. In §4, we show that the obtained equation reduces to the relativistic dissipative 
hydrodynamic equations in various frames with a choice of the macroscopic-frame 
vector, including the particle one and the energy one. Then, we compare our equa- 
tions with those proposed by Eckart, Landau-Lifshitz, and Stewart. In §5.1 and §5.2, 
we examine some properties of our equations, concerning transport coefficients and 
frames. In §5.3, we present a proof that all of our equations have the stable equi- 
librium state. The last section is devoted to a summary and concluding remarks. 
In Appendix A, we show that the most general form of the hydrodynamic equation 
as derived in the standard phenomenological way may admit the dissipative internal 
energy, pressure, and particle- number density. In Appendix B, after giving a brief 
account of the RG method using an example, we present a general ground of this 
method. In Appendix C, we present a detailed derivation of the first-excited modes 
in a generic local rest frame. 

§2. Preliminaries 

After a brief account of the basic properties of the relativistic Boltzmann equa- 
tion, we introduce a macroscopic-frame vector to specify a local rest frame on which 
the macroscopic dynamics is described: We shall see that the introduction of the 
macroscopic-frame vector enables us to have a coarse-grained and covariant equa- 
tion. As in the standard method like the Chapman-Enskog expansion method and 
others, 18 ^ we take the spatial inhomogeneity as the origin of the dissipation. We shall 
clarify the ad-hoc aspects of the ansatz's made in the standard Chapman-Enskog 
expansion and Maxwell-Grad moment methods for the derivation of the relativistic 
hydrodynamic equations. 

2.1. Relativistic Boltzmann equation 

As a simple example, we shall treat the classical relativistic system composed of 
identical particles *) . Then, the relativistic Boltzmann equation 18 ) for such a system 
reads 

P*d li f p (x) = C\f\ p (x), (2-1) 

where f p (x) denotes the one-particle distribution function defined in the phase space 
(x , p) with x^ being the space-time coordinate and p^ being the four momentum of 
the on-shell particle; p^ p^ = p 2 = m 2 and p° > 0. The right-hand side of Eq.(2-1) 
is called the collision integral, 

C[f] P (x) = — 4 S i ~o IP2, Ps) [f P2 ( x )fp 3 0*0 - f p (x)f Pl (x)^, 

pi 11 P2 ^ 2 P 3 3 

(2-2) 



An extension to multi-component systems is possible and will be presented elsewhere. 
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where uj(p , p\\p2 , P3) denotes the transition probability due to the microscopic two- 
particle interaction. To make explicit the correspondence to the general formulation 
given in Ref.33), we treat the momentum as a discrete variable, but the summation 
with respect to the momentum is interpreted as the integration; ^2 = fd^q with 
q being the spatial components of the four momentum q^. We remark that the 
transition probability u(p , p\\pi , Pz) contains the delta functions representing the 
energy-momentum conservation, 

w(p, P\\P2, Pz) oc (5 4 (p + pi -p 2 -P3), (2-3) 

and also has the symmetric properties due to the indistinguishability of the particles 
and the time reversal invariance of the microscopic transition probability, 

<*>(p, Pl\P2 , P3) = u{p 2 , P 3 \P, Pi) = W(P1 , P\P3, P2) = w(p 3 , P2\P1 , P)- (24) 

It should be stressed here that we have confined ourselves to the case in which the 
particle number is conserved in the collision process. 

The property of the transition probability shown in Eq.(2-4) leads to the follow- 
ing identity satisfied for an arbitrary vector (p p (x), 

£ ^ p {x)C[f] p {x) 4lE p lE^oE p TEj 

V %) V\ VI VZ 

XU>(p, Pl\p2,P3) (fp(x) + <Ppi(x) ~ <Pp 2 ( x ) ~ ^p 3 ( X )) 

X (/pa 0*0 fps ( x ) ~ fp( x ) fpi ( x )) ■ (2-5) 

A function (p p (x) is called a collision invariant (or summational invariant 18 ^) 
when it satisfies the following equation 

E 4^)C[/Ux) = 0. (2-6) 
p p 

As is easily confirmed by using the formula (2-5) and the property (2-3), (p p = 1 and 
p^ are collision invariants; 

E^[/1 P W = 0, (2-7) 

p 

^i/C[/] p (x) = 0, (2-8) 



p 

p 

which represent, of course, the conservation of the particle number, energy, and 
momentum by the collision process, respectively. We also see that the linear combi- 
nation of these collision invariants as given by 

i Pp {x) = a{x)+pi*b tx {x), (2-9) 

is also a collision invariant with a(x) and being arbitrary functions of x. This 

form is, in fact, known to be the most general form of a collision invariant. 18 ) 
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On account of Eq.'s (2-7) and (2-8), we have the continuity or balance equations 
for the particle current and the energy- momentum tensor T^ v ', 



d u T^{x) = d v 



v y 



= o, 



= 0, 



(2-10) 
(2-11) 



respectively. It is noted that while these equations have the same forms as the 
hydrodynamic equations, nothing about the dynamical properties is contained in 
these equations before the evolution of the distribution function f p (x) is obtained 
by solving Eq.(2-1). In the standard Chapman-Enskog expansion method, 18 ) these 
balance equations are rather used to obtain the time derivatives of the distribution 
function written in terms of the hydrodynamic quantities, order by order. 
The entropy current is defined by 



^(x)^-£^/ p (x)(ln/ p (x)-l). 



(2-12) 



Using the relativistic Boltzmann equation (2-1), the divergence of the entropy current 
reads 



(2-13) 



The above equation tells us that is conserved only if In f p (x) is a collision invariant, 
or a linear combination of the basic collision invariants (1 , p^) as 



In f p (x) = a(x) + p^ /3n(x), 



(2-14) 



with a(x) and /3 M (x) being arbitrary functions of x. In other words, the entropy- 
conserving distribution function is parametrized as 



f P (x) 



1 



(27T) 



cxp 



ji(x) -pV-u^x) 
T(x) 



(2-15) 



with u^{x) u fl (x) = 1, which is identified with the local equilibrium distribution 
function called the Juettner function 38 ) (the relativistic analog of the Maxwellian) : 
T(x), n(x), and u^(x) in Eq.(2-15) should be interpreted as the local temperature, 
the chemical potential, and the flow velocity, respectively. 

We note that for the local equilibrium distribution f p q (x) the collision integral 
identically vanishes, 



C[f\(x) = 0, 



(2-16) 



due to the energy-momentum conservation implemented in the transition probability 
(2-3). 
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2.2. Introduction of macroscopic- frame vector 

We are now in a position to introduce the key ingredient in the present work. 

Since we are interested in the hydrodynamic regime where the time and space 
dependence of the physical quantities are small, we try to solve Eq.(2-1) in the 
hydrodynamic regime where the space-time variation of f p {x) is small and the space- 
time scales are coarse-grained from those in the kinetic regime. To make a coarse 
graining with the Lorentz covariance being retained, we introduce a time-like Lorentz 
vector denoted by 

a**, (2-17) 

with 

a > 0. (2-18) 

Thus, a M specify the covariant but macroscopic coordinate system where the velocity 
field of the hydrodynamic flow is defined: Since such a coordinate system is called 
frame, we call the macroscopic-frame vector. Although a M may depend on the 
momentum p and the space-time coordinate x, i.e., 

a>+ = a£(x), (2-19) 

the time variation of it is supposed to be much smaller than that of the microscopic 
processes. We shall see that the separation of the scales between the kinetic and 
hydrodynamic regimes can be nicely achieved by the RG method. It should be 
stressed here that although a macroscopic vector is introduced also in the standard 
Chapman-Enskog expansion method, 18 ) the vector is identified as the flow velocity 
u^{x) from the outset. In our case, a,p(x) does not necessarily coincide with the flow 
velocity u^(x). 

Keeping in mind the above scale difference, we decompose the derivative <9 M into 
time-like and space-like ones in terms of the macroscopic-frame vector. Defining a 
projection operator to the space-like vector by 

a„(x) a v Ax) 

*?{x) = <r- M Q 2 (X ) (2-20) 



we have 



with 



glM _ a£{x)a u p {x) - a n (x )_?_ + JL (2-21) 



Yr^lix)^ X)d - ^ 
^- = A^{x)d v = Vr (2-23) 
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Then, the relativistic Boltzmann equation (2-1) in the new coordinate system 
(r , a^ 1 ) is written as 

p ■ a p (r , a) —f p (r , a) + p ■ Vf p (r , a) = C[f] p (r , a), (2-24) 

where a p (r , a) = a p (x) and f p {r , a) = f p (x). We remark the prefactor of the time 
derivative is a Lorentz scalar and positive definite; 

p ■ a p (r , a) > 0, (2-25) 

which is easily verified by taking the rest frame of p . 

Since we are interested in a hydrodynamic solution to Eq.(2-24) as mentioned 
above, we shall convert Eq.(2-24) into 

i-fp{r , a) = ) C[f] p (r , a) - e ) p ■ V/ P (r , a), (2-26) 

or p- a p (T , a) F p- a p (T , a) 

where a small quantity s has been introduced to express that the space derivatives 
are small for the system which we are interested in; e is called the non-uniformity 
parameter 18 ) and may be identified with the ratio of the average particle distance 
over the mean free path, i.e., the Knudsen number. 

Since e appears in front of the second term of the right-hand side of Eq.(2-26), 
the relativistic Boltzmann equation has a form to which the perturbative expansion 
is applicable. In fact, this form of the Boltzmann equation but with a p = is also 
the starting point for the standard Chapman-Enskog expansion method, 18 - ) together 
with an ad-hoc ansatz on the order of the time derivative of the distribution function. 
We stress that this seemingly mere rewrite of the equation has a physical significance; 
it expresses a natural assumption that only the spatial inhomogeneity over distances 
of the order of the mean free path is the origin of the dissipation. We remark that the 
RG method applied to non-relativistic Boltzmann equation with the corresponding 
assumption leads to the Navier-Stokes equation; 25 )' 27 )' 34 )' 35 ) the present approach is 
simply a covariantization of the non-relativistic case. 

For setting up the perturbative expansion in a consistent way with the physical 
picture of the origin of the dissipation, we shall take the coordinate system where 
a p {r , a) has no r dependence, i.e., 

a£(T,a) = a£(<7). (2-27) 

Then, Eq.(2-26) takes the following form, 

4-fpir , a) = l — C[f] p (T ,a)-e ^ p • V f p (r , a). (2-28) 

<9r p- a p {a) F p ■ a p (a) 

Now let q p (cr) is a physical quantity of a particle with a momentum p M at a 
position cj m ; then the total amount of the quantity is given by 

Q P (r) = Jd 3 aq p (a)f p (T, a). (2-29) 
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When q p (a) = p • a p (a), the time variation of Q p (t) is given by 

-^Q p {t) = d 3 ap- a p (a) —f p (T , a), 

= -ep» jd 3 a^Lf p (r, a) + jd 3 a C[/] p (r , a), 

= jd 3 aC[f] p (T,a), (2-30) 

where we have used Eq.(2-28) and neglected the contribution from the spatial bound- 
ary of the system. Here, we consider a trajectory of one particle in the interval 
between each collisions: Substituting C[f] p (r , a) = into Eq.(2-30), we find that 
Q p (t) is a conserved quantity and q p (o~) is a corresponding density. Thus, the form 
of the relativistic Boltzmann equation as given in Eq.(2-28) tells us that the physical 
quantity that is transported by each particle in the system is given by 

p-a p {a). (2-31) 

It is to be noted that we can control what is the flow represented in our theory 
by varying the specific expression of a p (a). Because of this freedom inherent in 
our coordinate system, our theory may lead to various hydrodynamic equations, 
including the ones in the energy and particle frames for non- ideal fluids. 

2.3. A brief description of the standard methods 

Before developing our analysis based on the RG method, we briefly summarize 
the ad-hoc aspects in the standard methods such as the Chapman-Enskog expansion 
and Maxwell-Grad moment methods. 

In the standard Chapman-Enskog expansion method, one starts from the fol- 
lowing form of the Boltzmann equation, 

V ■ u DcEfp(x) = -ep- VcEfpix) + C[f} p (x), (2-32) 

where D C e = and V£ E = Ag E d v with Ag E = - u>* u v . Then, one makes 

the perturbative expansion 

f p (x) = fj°\x) +eff\x) + e 2 ff\x) + ■ ■ ■ = /f(x) + 5f p (x), (2-33) 
DceU(x) =eD C Ef ( p l \x) + e 2 Dce^Hx) + ■■■ . (2-34) 

Then, it is found that the zeroth-order solution is given by the local equilibrium 
distribution function, i.e., the Juettner function given by Eq.(2-15); 

f^(x) = f^(x). (2-35) 

It is customary, 18 ) without any physical foundation, to assume that the particle- 
number density and internal energy in the non-equilibrium state is the same as those 
in the local equilibrium state, and set as follows; 



v ' ' 



(2-36) 
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e = u. 



E^4 0) 



For consistency, one also imposes the constraints to the higher-order terms 



8n = u L 



5e = u 



o. 



u„ = 0. 



(2-37) 



(2-38) 



(2-39) 



To obtain the hydrodynamic equation in the particle (Eckart) frame, another con- 
straint is imposed 



CE 



E -pPvtu 



0. 



(2-40) 



Instead, if one wants to obtain the hydrodynamic equation in the energy (Landau- 
Lifshitz) frame, one also imposes the constraint 



^CE ^ vp 



4& 



E ^p»pp 6 fp 



0. 



(2-41) 



The constraints imposed to the distribution function in the higher orders are called 
the conditions of fit; the zeroth-order constraint is not a constraint but an identity. 

It is noteworthy that although a foundation for them has never been given, 
such conditions of fit (2-36)-(2-41) are also imposed in an ad-hoc way even when 
the Maxwell-Grad moment method is adopted. 18 ) We stress here that these con- 
straints are actually equivalent with a strong physical assumption that there are no 
particle-number density nor internal energy of the dissipative origin, although the 
distribution function in the non-equilibrium state is quite different from that in the 
local equilibrium state. It is, therefore, an urgent but yet unsolved problem to verify 
or elaborate these ad-hoc constraints somehow, say, from a microscopic theory, or 
by experiment, if possible. 

In the RG method which we adopt, one needs no such conditions of fit for deriva- 
tion, and rather the correct forms of them are obtained as a property of the derived 
equation once the frame is specified by the macroscopic-frame vector: We shall see 
that the conditions of fit in the energy frame is compatible with the underlying 
Boltzmann equation and physical, but those in the particle frame is not and thus 
will never be satisfied in any physical system. 



§3. Reduction of Relativistic Boltzmann Equation with 
Renormalization-group Method 

In this section, starting from Eq.(2-28) we shall derive the relativistic dissipative 
hydrodynamic equation as the infrared asymptotic dynamics of the classical rela- 
tivistic Boltzmann equation on the basis of the RG method: The five hydrodynamic 
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variables, i.e., the flow velocity, local temperature, and particle-number density (or 
chemical potential), naturally correspond to the zero modes of the linearized colli- 
sion operator. Then, without recourse to any ansatz such as the conditions of fit 
given by Eq.'s (2-36)-(2-41), the excited modes which are to modify the local equilib- 
rium distribution function are naturally defined in the sense that they are precisely 
orthogonal to the zero modes with a properly defined inner product for the distribu- 
tion functions. It will be shown on the basis of the inner product that the so-called 
Burnett term does not affect the hydrodynamic equation owing to the fact that the 
hydrodynamic modes are the zero modes of the linearized collision operator. 

3.1. Construction of the approximate solution around arbitrary initial time 

In accordance with the general formulation of the RG method 29 )' 30 -*' 32 )' 33 ) *\ 
we first try to obtain the perturbative solution f p to Eq.(2-28) around the arbitrary 
initial time r = To with the initial value / p (to, cr); 

f P {r = r , a; r ) = / p (r , a), (3-1) 

where we have made explicit that the solution has the To dependence. The initial 
value is not yet specified, we suppose that the initial value is on an exact solution. 
The initial value as well as the solution is expanded with respect to e as follows; 

f P (r ,a;r ) = ff^r ,a;r ) + e ff\r , a ; r ) + e 2 # 2 >(t , a ; r ) + • • • , 

(3-2) 

and 

f p (r , a) = 4°)(r , a) + £ /«(t , a) + e 2 f\r , a) + • • • . (3-3) 
The respective initial conditions at r = To are set up as 

(t , a; T ) = ff (t , a) for I = 0, 1, 2, • • • . (3-4) 

In the expansion, the zeroth-order value fp°\ro , a ; tq) = fp \ro , cr) is supposed to 
be as close as possible to an exact solution. 

Substituting the above expansions into Eq.(2-28) in the T-independent but To- 
dependent coordinate system with 

o£(t, a) = a^(T , a) = a$(a; t ), (3-5) 

we obtain the series of the perturbative equations with respect to e. 
Now the zeroth-order equation reads 

(t, a; to) = } -C[f^] p (r,a;r ). (3-6) 

or p p- a p (cr ; r ) 

Since we are interested in the slow motion which would be realized asymptotically 
as t — y oo, we should take the following stationary solution or the fixed point, 

^/f(T,a;T ) = 0, (3-7) 
See Appendix B for a brief but self-contained account of the RG method. 
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which is realized when 



P ■ a p {cr ; r ) 



C[f\(r,a;T ) = 0, 



(3-8) 



for arbitrary a. Thus, we see that ln/^ 0) (r , a; tq) can be represented linear 
combination of the five collision invariants (1 , p^) as mentioned in the last section, 
and hence fp°\r , a ; To) is found to be a local equilibrium distribution function and 
thus given by the Juettner function (2-15): 



(27T) 



cxp 



MgJ t ) -p^u^a; t ) 
T(a; T ) 



= / p e >;r ). (3-9) 



with u M (<7 ; to) u^a ; to) = 1, which implies that 

4°)(t , a) = /f (t , a ; tq) = / p c > ; t ). 



(3-10) 



It should be noticed that the five would-be integration constants T(a ; To), [i(cr ; To), 
and u^(a; To) are independent of r but may depend on To as well as a. In the 
following, we shall suppress the coordinate arguments (a ; to) and the momentum 
subscript, e.g., p when no misunderstanding is expected. 

3.2. The linearized collision operator and inner product 
Now the first-order equation reads 

d_ 



-f p [1) (T) = ^A p JV(T)+F p , 



(3-11) 



where the linear evolution operator A and the inhomogeneous term F are defined by 

1 d 



and 



Apq — rsr C[/]p 



P ■ a p df q 

= ~ 777 ElUlI! ^o ^(p, pi|p 2 , P3) 

p ■ a v 2! pi ^ rfi ^ b, 

A \ u P2qJp 3 ' Jp2 °P3<? °pqJ Pl Jp v piq)i 



F P = ~ 



p-a p 



(3-12) 



(3-13) 



respectively. 

To obtain the solution which describes a slow motion, it is convenient to first 
analyze the spectral properties of A. For this purpose, we convert A to another 
linear operator 



L = f^ 1 A / cq , 



(3-14) 
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with the diagonal matrix fpq = fp q 5 pq ; the explicit form of L is given by 

Lpq = ^~7T 91 E Tfl E ~fi E ~fl ' Pl \ P2 ' **) ( S P* + S P19 ~ S P2Q ~ <W 

P-a p z. pi Pl p2 p 2 ^ p 3 

(3-15) 

Here, we have used the identity 

"(p , Pi\P2 , p 3 ) / P ° 2 q / P eq = w(p , Pi|P2 , Pa) / P eq /£, (3-16) 

which follows from Eq.'s (2-3) and (3-9). 

Let us define the inner product between arbitrary non-zero vectors (p and ip by 

V) = E ^(p-a p )fp<p p ^ p . (3-17) 
p P 

We note that the norm defined through this inner product is positive definite 

(^^)=^i(p-a p )/;fe) 2 >0 for <p p ^0, (3-18) 
p P 



since 



p-a p >0 (3-19) 

in accord with Eq.(2-25). Notice that the other factors than p ■ a p in Eq.(3-18) are 
all positive definite. We shall see that this positive definiteness (3-18) of the inner 
product plays an essential role in making the resultant hydrodynamic equations 
assure the stability of the thermal equilibrium state, as it should be, in contrast to 
some phenomenological equations. 

With the inner product, it is found that L is self- adjoint 

V Pl 1 P2 Z P3 6 

W(P , Pl\P2 , P3) /pi 1 (fp + f Pl ~ fp 2 ~ Ppz) (YV + Ippi - ^ P 2 ~ V'ps) 

= (L^,V), (3-20) 
and semi-negative definite 

LLP) = "2! 4 E ^ E E E i 

P Pl 1 P2 Z P3 6 

xuj( p , pilpa , p 3 ) f° q f£ (f p + if Pl - (f P2 - f P3 ) 2 
< 0, (3-21) 

which means that the eigen values of L are zero or negative. In the derivation of 
Eq's (3-20) and (3-21), we have used Eq.'s (2-3) and (3-16). 
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The eigen vectors belonging to the zero eigen value are found to be 

( jf for a = u, 

| lxm for a = 4, 

which span the kernel of L and satisfy 

[L^] p = 0. (3-23) 

We call <Pq the zero modes. It is noted that these zero modes described by the five 
vectors are collision invariants shown in Eq.'s (2-7) and (2-8), and the factor m in 
ip\ v is introduced merely for convenience so that our method can be applied to the 
case of massless particles. 

Following Ref.33), we define the projection operator Pq onto the kernel of L 
which is called the Po space and the projection operator Qq onto the Qo space 
complement to the Po space: 

[Pi>^] p = 'PdpV^(^,^), (3-24) 
Qo = l- Po, (3-25) 

where r/jj^g is the inverse matrix of the the Po-space metric matrix r^f defined by 

r^ = (^ lV g). (3-26) 

3.3. First-order solution 

The solution to Eq.(3Tl) with the initial condition /^(t = To) = i.e., 
/p^(r = r , a ; r ) = fp\(J ; r ) is expressed as 

/(l) (r) = e (r~r )A + A -l Q o p j + (r _ Tq) p o ^ _ A -l g Q ^ (3.^ 

where we have introduced the modified projection operators 

P = / eq Po/ eq -\ (3-28) 
Qo^nQon' 1 . (3-29) 

We remark that the first term in Eq.(3-27) would be a fast motion coming from the 
Qo space, which can be simply eliminated by choosing the initial value /W, which 
has not yet been specified, as 

f W = P ) (T ) = -A~ 1 Q F. (3-30) 

Thus, we have the first-order solution 

/(i) ( T ) = {j- T ) Po F - A" 1 Q F, (3-31) 

with the initial value (3-30). We notice the appearance of the secular term propor- 
tional to t — tq, which apparently invalidates the perturbative solution when \t — tq| 
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becomes large. It is worth mentioning that the standard Chapman-Enskog expansion 
method includes a set of conditions for not making secular terms appear; the condi- 
tions are the solvability conditions of the balance equations (2-10) and (2-11). 18 ) In 
applying the solubility conditions, one needs apply the ad-hoc constraints on the dis- 
tribution function for defining the flow (eg. Eckart flow or Landau-Lifshitz flow) as 
well as the ad-hoc constraints on the particle-number density and internal energy, 18 ) 
as given by Eq.'s(2-36)-(2-41). In the present RG method, secular terms are allowed 
to appear and no constraints are imposed on the distribution function; rather the 
secular terms will be utilized to obtain the slow dynamics. 

We remark here that we could apply the RG method here to Eq.(3-31), which will 
give the relativistic Euler equation without dissipation effects. To get a dissipative 
hydrodynamic equation, we need to proceed to the second order in our method. 

3.4. Second-order solution 

The second-order equation is written as 

to = E A pi fV to + (t-t ) 2 G p + (t- to) H p + I p , (3-32) 

q 

where 



Bpqr ~ p-a p 2df q dfr mp 



f=f e 



1 1 1 



J2 i E i 12 tp "(p , Pl \ P2 , P3 ) 



p ■ a v 2 2! ^ p\ ^ p% ^ p% 

p pi f-r P2 fj. P3 i.i 

X ($p2q 0~p3r ~i~ &p2T 3psq ®~pq 0~pir 0~pr ^piq)i (3 - 33) 

G P ^EE B ™r [Po F] q [Po F] r , (3-34) 

q r 

ff^-EE *W ( [Po F] q [A- 1 Q F] r + [A' 1 Qo F] q [P F] r ) 

q r 

-rV^ V ^^] p , (3-35) 

jj ■ Lip 

^EE V [A' 1 Qo F] [A- 1 Q F] r + — p-W[A- l Q Q F] (3-36) 

q r 

The solution to Eq.(3-32) is found to be 

/(2)( T ) = e ( T - T ^ A [f^ + 2 A- 3 Qo G + A" 2 Q H + A' 1 Q /} 
+ i (r - r ) 3 Po G + \ (r - r ) 2 {P H - 2 A' 1 Q g] 
+ (r - T ) {Po I - 2 A~ 2 Qo G - A- 1 Q p} 
+ { - 2A- 3 QoG-A- 2 Q H -A^Qol}. (3-37) 
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Again the would-be fast motion can be eliminated by a choice of the initial value 
/( 2 ) so that 

/( 2 ) = /( 2 ) (r ) = -2 A~ 3 Q G - A' 2 Q H- A' 1 Q 1. (3-38) 
Then, we have a second-order solution 

/ {2) (r) = ^ (r - to) 3 Po G + ^ (r - to) 2 {p ^ - 2 A" 1 Qo g} 
+ (r - r ) {Po / - 2 Qo G - A' 1 Q p} 
+ { -2A- 3 Q G-A- 2 Q H - A^Qol}, (3-39) 

We notice again the appearance of secular terms and that no constraints on the 
solution are imposed for defining the flow, in contrast to the standard Chapman- 
Enskog expansion method. 18 ) 

Summing up the perturbative solutions up to the second order, we have an 
approximate solution around r ~ To to this order; 

fp(r , a ; to) = ff\r , a ; r ) + e ff\r , a ; r ) + e 2 ff\r , a ; r ) + 0(e 3 ). 

(3-40) 

We emphasize that this solution contains the secular terms which apparently invali- 
dates the perturbative expansion for r away from the initial time To- 

3.5. RG improvement of perturbative expansion 

The point of the RG method lies in the fact that we can utilize the secular terms 
to obtain an asymptotic solution valid in a global domain. Now we may see that 
we have a family of curves f p {r , a ; To) parameterized with To- They are all on the 
exact solution f p (a; r) at r = tq up to 0(e 3 ), but only valid locally for r near To- 
So it is conceivable that the envelope of the family of curves which contacts with 
each local solution at r = To will give a global solution in our asymptotic situation. 
According to the classical theory of envelopes, the envelope which contact with any 
curve in the family at t = tq is obtained by *) 



^-; p (t,^;t ; 



= 0, (3-41) 

T()=T 



or explicitly 
d 



— {r-sA- 1 Q F}-eP F 



-e 2 {p I-2A~ 2 Q G-A- 1 Q H}+O(e 3 ) = 0. 



(3-42) 



See Appendix B for the foundation of this procedure. 
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This envelope equation is the basic equation in the RG method and gives the equation 
of motion governing the dynamics of the five slow variables T(a; r), /j,(a ; r) and 
u^(a ; r) in fp q (a ; r). The global solution in the asymptotic region is given as an 
envelope function, 

/e p (t, cr) = f p {r, ct;to = t) 

= f^(a-T)-e[A~ 1 Q F] p (a;r) 

-e 2 {2[A- ;i QoG} p (a;T) 

+ [A' 2 Q H] p (a; r) + [AT 1 Q l] p (a ; r)} + 0(e 3 ), (3-43) 

where the exact solution of Eq.(3-42) is inserted. As is proved in Appendix B, 
the envelope function /ep(t , a) satisfies Eq.(2-28) in a global domain up to 0(e 3 ) 
owing to the condition (3-41), although f p {r , a ; To) itself was constructed as a local 
solution around r ~ To- Thus, one sees that /e p (t , cr) now describes a coarse- 
grained evolution of the one-particle distribution function in Eq.(2-28), because the 
time-derivatives of the quantities in /e p (t ; a ) are au i n the order of e or higher. We 
emphasize that we have derived the slow-motion equation of Eq.(2-28) in the form 
of the pair of Eq.'s (3-42) and (3-43). 

3.6. Reduction of the RG equation to generic hydrodynamic equation 

Now let us see that the RG/Envelope equation (3-42) is actually the hydrody- 
namic equation governing the five slow variables, T(cr; r), (i(a ; r), and u^(a; r). 
To show this explicitly, we apply Po from the left and then take the inner product 
with the five zero modes ify. In this procedure, we first note that the direct use of 
the definitions of Eq.'s (3-24), (3-26), and (3-28) leads to the following identity; 

E i (p ■ %>) ^p [^4 = Ei(f a p) ( 344 ) 

p p 

Furthermore, noting that y?g are the collision invariants as shown in Eq.'s (2-7) and 
(2-8), we have 

E i (P ■ ^ E E B w <5o F] g [A" 1 Qo F] r = 0. (3-45) 

p q r 

Here, we have used the relation 

E B Pir [A' 1 Qo F] q [A" 1 QoF] r = -L- C[A- 1 Q F] p , (3-46) 

q r P p 

which follows from the definitions of Eq.'s (2-2) and (3-33). 
Thus, we have 

j/eq _ £ [ A -l Q Q F Q + 0(£ 3) = 0, (3-47) 



E^o P 



(p-a p ) — + ep-V 
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Putting back e = 1, we arrive at 



« = o 



(3-48) 



with 



J lst - 2^ „0 P 



{/ P cq -[^ 1( 3oF] p }, 



(3-49) 



p 



where we have used 



(p-a p ) — +p- V=p fM 8^. 



(3-50) 



It is noted that Jf s " perfectly agrees with the one obtained by inserting the solution 
hp(r, a) in Eq.(3-43) into W and in Eq.'s (2-10) and (2-11): 



Therefore, we conclude that Eq.(3-48) is identically the relativistic dissipative hy- 
drodynamic equation: We note that the subscript of "1st" in Jf s " means that the 
obtained equation is in a class of the so-called first-order relativistic dissipative hy- 
drodynamics. As is shown in Appendix A, the first-order equations can be derived 
phenomenologically with use of the entropy current which includes dissipative effects 
up to the first order. 

It is noted that the hydrodynamic equation (3-48) with the currents (3-49) still 
contains the macroscopic- frame vector dp, which is now dependent on r as well as on 
a. We shall see that Eq.(3-49) reduces to the currents in various frames with a choice 
of ap 1 . In other words, we have obtained the relativistic hydrodynamic equations for 
a viscous fluid in a generic frame, which is a kind of the master equation from which 
various relativistic hydrodynamic equations are deduced. As far as we are aware of, 
this is the first time when such a generic hydrodynamic equation for the relativistic 
viscous system is obtained. This is one of the main results in the present work. We 
stress that this was made possible because any ad-hoc ansatz such as the conditions 
of fit is not necessary in the RG method. 

It is also worth mentioning that the problematic Burnett term is absent in 
Eq.(3-47) thanks to Eq.(3-45). If the Burnett term were to remain, the particle- 
number and energy-momentum conservation laws are lost; moreover, boundary con- 
ditions might have to be taken care of simultaneously because its magnitude is com- 
parable to the Burnett term. 39 ) In fact, the presence of the Burnett term is known 
to be inevitable when the Chapman-Enskog expansion method 39 ^ is applied to derive 
the Navier-Stokes equation from the non-relativistic Boltzmann equation *). 

*-* We did not recognize this important point in Ref.36) that the RG method naturally gives 
a hydrodynamic equation free from the Burnett term. We expect that the mechanism similar to 
Eq.(3-45) is realized in the non-relativistic case. 




(3-51) 
(3-52) 
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We decompose Jf s " into two parts as Jf s " = j[gl^ a + SJ^., where 

V 

Kst = ~ E 75? ^ [A' 1 QoF] =-(#?, L- 1 Q /i-i F >, (3-54) 



p 



with 



(3-55) 

P Up 

Needless to say, j[^ a and 5Jf s " represent the currents in the perfect-fluid and 



dissipative part, respectively. Corresponding to J^l 1 "* an d ^lsti ^ M an d 1 
Eq.'s (3-51) and (3-52) are decomposed as 

N» = NW» + <W*, (3-56) 

T^ = T^ U + 5T^, (3-57) 

where 

N^^m- 1 ^, (3-58) 

Mf^m -1 ^, (3-59) 

T(°)^ = J^"", (3-60) 

ST^ = SJ^f (3-61) 

For later convenience, we present a simpler form of 5J^. First, we note that F 
in Eq.(3-13) is expressed as 

F P = -f^{K m ~ 1 ^4-^ v ^) = - f ^^^ (3 ' 62) 

where 

f -V„(u„/T) for a = i/, 
= < , 3-63 
^ \ m' l V^/T) for a = 4. 

Using the above representation of F and the identity 

{i P ,L- 1 Q ^) = {Q ip,L- 1 Q iP), (3-64) 
we can reduce the dissipative part SJ^ to the following form: 

6Jg = vr*X vP , (3-65) 

where 

< a ^ = (^ a ,L-Vf)- (3-66) 
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Here, we have introduced an important new vector defined by 



Pip 



(3-67) 



We call ip^p the first- excited modes. It is noteworthy that <5Jf s " is represented as a 



product of rji aul3 and X v q: r(^ vp has some information about the transport coeffi- 
cients, while X v p is identical to the corresponding thermodynamic forces. 

Finally, we write down the relativistic dissipative hydrodynamic equation defined 
in the covariant coordinate system {a ; r) : 



jiavfi 



v d 



/p -(l + [L-Vf] p ^) 



0, (3-68) 



which can be derived straightforwardly from Eq.'s (3-47) and (3-62). This form of 
the relativistic dissipative hydrodynamic equation will be found to play an essential 
role in the stability analysis to be presented in §5.3. 

§4. Relativistic Hydrodynamic Equation with Specific 
Macroscopic-frame Vectors 



In this section, we give explicit forms of the relativistic dissipative hydrodynamic 
equations by integrating out the right-hand side of Eq.'s (3-53) and (3-54) with 
respect to the momentum p^. 

A remark is in order here: We calculate and present the thermodynamic quan- 
tities and transport coefficients using our model equation, i.e., the relativistic Boltz- 
mann equation, for completeness. Then, the explicit forms of them are inherently 
for the relativistic rarefied gas. We would like to remind the reader, however, that 
the main purpose of the present work is to determine the form of the relativistic 
hydrodynamic equations for a viscous fluid, and expect that the forms of the macro- 
scopic hydrodynamic equations which contain the thermodynamic quantities and 
transport coefficients only parametrically, and hence the forms are independent of 
the microscopic expressions of these quantities. 



After the integration, the currents of perfect-fluid part j[^ a 

. pA i»> =r (0)F for a = I/j 

: m 7V (0)/i for a = 4, 



reads 



J 



(0)/ia 
1st 



ev?u v 



(4-1) 



mnw 



where e, p, and n denote the internal energy, the pressure, and the particle-number 
density, respectively. These quantities are defined by 



n = Y 

p 



p 



(27T) 



4?rm er z 



f; q (p • uf 



ran 



K 2 (z) 



— z 



^(-l/3p"/V)=nT, 



(4-2) 



(4-3) 



(4-4) 
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where we have introduced the second- and third-order modified Bessel functions 
K2(z) and K^(z) with z being the dimensionless variable defined by 

(4-5) 

The explicit form of the modified Bessel functions is presented in Eq.(C-8) in Ap- 
pendix C. 

4.1. A generic choice of the macroscopic- frame vector 

The dissipative part of the currents <5Jf s " depends on the macroscopic-frame 
vector dp explicitly, in the contrast of J^l*"*- Here, we shall consider <5Jf s " with a 
choice of ap 1 . As a simple but nontrivial choice, let us take the following set of the 
macroscopic- frame vectors with 9 being a constant; 

a£ = — ( (p • u) cos 9 + m sin d) = 6%. (4-6) 
v p- u V / p 

In Eq.(4-6), we note that the factor m is introduced simply to make the expression 
dimensionless, so our method is also applicable to the case of massless particles. 
A remark is in order here: The inequality (2-18) of a° leads to a restriction on 9 

when ap 1 = Op in Eq.(3T9), as (p-u) cos9 + m sin 9 = y 7 (p ■ u) 2 + in 2 sm(6 + x) > 
with tanx = (p ■ u)/m > 1, which implies that 

7T „ IT 

-j<0<-, (4-7) 

because 7r/4 < x < 7r/2. 

Under these settings, we shall now give the explicit representation of <5J{^ = 
^ixavp x v p_ First, we consider in Eq.(3-65). Inserting the equality ap 1 = 9p, we 
have 

A»» = g»» - u »u u = A» u , (4-8) 
so that the differential operator with respect to <r M reads 

Notice that A^ v and V M are now identical to the familiar projection matrix onto 
the subspace complement to u M and the covariant spatial differential operator, 1S ) 
respectively Thus, the thermodynamic forces X^ a are reduced to 

( -VJu v /T) for a = v, 
I™ = <^ , 4-10 
M \ m- 1 V M (/x/T) for a = 4. 

Next, we derive an explicit form of r/^ ai//3 in Eq.(3-65), which task is tantamount 
to obtaining an explicit representation of ip^ . After a straightforward calculation 
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presented in Appendix C, one can find that ip^^ are written as 

+ J£ Y 3 (9) u u + r p Y 3 {9)u» + ft? for a = v, 
, fl p Z x (9) u» + J£ Z 2 {9) for q = 4. 



ua 
<Plp 



Here, we have introduced lip, Jp, and defined by 



(Hp, Tf^) = (J7 P , J£, tt: 



(4-11) 



(4-12) 



and Yi(0 
Y 1 (9 

Y 2 (9 

Y 3 (9 
Zi{9 

Z 2 {9 



n p={\~ 7) (P • + ((7 - 1) T h - 7 r) (p • u) - \ m 2 , (4-13) 
J^ = -({p-u)-Th)A^p v , (4-14) 
< = ^> pP(T , (4-15) 

, Y 2 {9), Y 3 (9), Z x {9), and Z 2 {9) defined by 

-3z 2 sin 2 (9 

_ z 2 cos 2 + z 2 (37-4) sin 2 0-3z[l-(/i-l)(7-l)] cos#sin#' 



(4-16) 



z 2 cos 2 9 — z 2 sin 2 1 



z 2 cos 2 6> + z 2 (37-4) sin 2 6»-3z[l-(/i-l)(7-l)] cos 9 sin 6* 



-z sint 



h cos + 2 sin 9 



(4-17) 
(4-18) 



3z 2 sin# cos# 



z 2 cos 2 9 + z 2 (37-4) sin 2 0-3z[l-(/i-l)(7-l)] cos 9 sin 9 



z cos 9 



h cos 9 + z sin 



(4-19) 
(4-20) 



respectively. It is noted that IJ p , Jp, and ~Kp V belong to the Qo space. n p , Jp, and 
itp U are independent of 9, which are identically the microscopic representations of 
dissipative currents in the literature. 18 ) The definitions of h, 7, and A^f" 7 used in 
Eq.'s (4-13)-(4-20) are 



h = 



e + p 



nT ' 

7 = 1 + (z 2 - h 2 + 5h - I) -1 , 

Alivpa = \( AjJ-P AM 1 /W /\J/p _ ^ . MI/ /ipoA 
2 V 3 /' 



(4-21) 
(4-22) 
(4-23) 
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where h and 7 denote the reduced enthalpy per particle and the ratio of the heat 
capacities, respectively. 

With use of ^ip in Eq.(4-ll), we can write down r] paul3 = { <^ Q , L^ 1 tpf ) as 

rf x pV(T = -T( (y^O) u p u p - Y 2 (0) A pp ^j (Yi(0) u v u u - Y 2 {9) A uu ^j 

+ T 2 \ Yf(6) (u p u v A pa + u p u a A pv + u p u v A pa + u p u a A pv ) 
-2Tr]A ppua , (4-24) 

rft p u4 = rffw = -T( (y^O) u p u p - Y 2 {6) A pp ^j Z x (6) u v 

+ T 2 \Y 3 (0) Z 2 (6) (u p A pv + u p A pu ), (4-25) 
rtf 4 " 4 = -TC,Z\(d)u p u v + T 2 \Z 2 {6)A pu . (4-26) 

Here, we have introduced the transport coefficients, i.e., the bulk viscosity £, the 
heat conductivity A, and the shear viscosity 77 given by 

= — L^H), (4-27) 

A^^(J>,L- 1 J M ), (4-28) 

V^-j^i^^L- 1 ^), (4-29) 

respectively. In §5.1, we show that these quantities are independent of 0. In the 
derivations of Eq.'s (4-24)-(4-26), we have used the following identities 

( > , L- 1 J v ) = \{ J a , L- 1 J a > A pu , (4-30) 

( 7f^ , L- 1 TT pcr ) = — ( n ab , L- 1 n ab ) A^ vpa . (4-31) 
5 

Now let us rewrite the expressions of £, A, and r\ in a more familiar form, i.e., the 
Green-Kubo formula 40 ) in the linear response theory. 24 )> 41 ) With use of the identity 

poo 

E i L ~\ g o = - / ds E l eS % g o> ( 4 - 32 ) 

q J ° q 

we can rewrite Eq's (4- 27)- (4- 29) as 

poo 

C= ds R c (s), (4-33) 
J 

o 

ds R x (s), (4-34) 

poo 

v = ds R v (s), (4-35) 



where 

Rs(s) = 

T 



R c (s) = ^(n(0),II(s)}, (4-36) 
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Rx(s) = -^ (J>(0), J^s)), (4-37) 

^W=1^<^(0),^W>, (4-38) 

with 

(n p (s), jjt{8), = 2 t^U J « > O- ( 4 " 3g ) 

<? 

It is noted that R$(s), R\(s), and R v (s) in Eq.'s (4-36)-(4-38) are called the relaxation 
function in the linear response theory. 24 )> 41 ) 

Finally, the dissipative currents SJ^ are obtained from r^ av ^ in Eq.'s (4-24)- 
(4-26) and X uf3 in Eq.(4-10), as follows: 



0J lst 



- T A y 3 (0) («" X} + u" X^j) + 2 r? = 5T» V for a = v, 
QZ x (Q)v?X n -T\Z 2 {0)X» = m5N^ for a = 4, 

(4-40) 

where the new thermodynamic forces Xn-, Xj, and are defined by 

X n = -T (Yi{9) u" u v - Y 2 (6) A^ X pv - T Z X {Q) u» X p4 

= -Y 2 (9)V-u, (4-41) 

X » = -T A^ (Y 3 (6) u"(X vp + X pv ) + Z 2 (9) X„ 4 ) 

= r 3 (9)TV,(l/T) - Z 2 (8) z- 1 V M (/i/T), (4-42) 



= -TA^ pu X pa 



(4-43) 



Here, the following relations have been used: X pa = and u v V p u u = 0. Thus, 

we have arrived at the generic form of the currents = j[sl^ a + ^^ist f° r an 
arbitrary 6 specifying the local rest frame of the flow. Equations (4-1) and (4-40) 
together with Eq.'s (4-41)-(4-43) are the main results in this section. 

To have a physical intuition into the effects of the dissipation on the hydrody- 
namic equations, we also decompose 5Jf s " in Eq.(4-40) into various tensors 

{ SevP u v - 8p Ai" + Qi* u v + Of vP + III™ = 5T pu for a = u, 
6J& = \ (4-44) 
] m5nu p + mv 11 = m5N^, for a = 4, 



where 



Se = 5JZ u p u v = 5T» V u p u v , (4-45) 
Sp = -1/3 5JZ A pv = -1/3 ST"" zV, (^) 
Q» = SJZ UvAf = 5T"> u v Af, (4-47) 
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= 5J&. A p f v = ST pa Apf, (4-48) 
5n = m _1 d J{f t u p = 5N^ Uli , (4-49) 
v» = m^SJ^ A/ = 5N U A/. (4-50) 

Here, Se, Sp, and 5n represent the contribution of the dissipation to the internal 
energy, pressure and particle-number density, respectively. We stress that the dissi- 
pation in the relativistic system generically gives rise to an additional contribution 
to such would-be thermodynamic quantities, as well as the dissipative currents, Q^, 
II^ U , and v^. Furthermore, it is to be noted that the existence of be defined in 
Eq.(4-45) is directly related with the definition of the local rest frame as is seen from 
Eq.(l-l). 

Inserting Eq.(4-40) into Eq.'s (4-45)- (4- 50), we have 
5e = (Y 1 (9)X n 

= -(Y 1 (6)Y 2 (e)V-u, (4-51) 
5p = tY 2 {9)X n 

= -CY 2 2 (6)V-u, (4-52) 

= -T\Y 3 {e)xy 

= —TX (^(ejTV^l/T) - Y 3 (0) Z 2 {6) z- 1 VOi/T)) , (4-53) 
= 2rjX% u 

= 2 v A^f" T V p u a , (4-54) 
m5n = (Z 1 (6)X n 

= -(Y 2 (6)Z 1 (e)V-u, (4-55) 
mv» = -T\Z 2 {6)X^ 

= —TX (Y 3 (0) Z 2 {9)TV^{l/T) - Z%(0) z' 1 VV/T)) . (4-56) 

Equations (4-51)-(4-56) also constitute main results in this section, and we call them 
the constitutive equations. 

We shall consider here the constitutive equations only for a few values of 6 
which satisfies Eq.(4-7), but show that our generic form of the constitutive equations 
contains ones in the two well-known frames, i.e., the energy one and the particle one, 
which were introduced in §1. 

4.2. Energy frame 

With the simplest choice 6 = 0, i.e., dp = we have 

be = 0, (4-57) 

5p=-(V-u, (4-58) 

= 0, (4-59) 

= 2 r]A^P a V pUcT , (4-60) 

mSn = 0, (4-61) 
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m / = A™V^, (4-62) 

where we have used the following relations obtained from Eq.'s (4-16)-(4-20), 

n(0) = 0, (4-63) 

F 2 (0) = 1, (4-64) 

F 3 (0) = 0, (4-65) 

Zi(0) = 0, (4-66) 

Z 2 (0) = z/h. (4-67) 

It is noted that the choice 9 = gives Q^ 1 = as well as / 0, which tells us that 
our hydrodynamic equation with = is the one in the energy frame. We remark 
that the dissipation gives rise to an additional pressure 5p while the dissipative 
internal energy Se and particle-number density 5n are absent in this frame. Thus, 
the conditions of fit (2-38) used in the usual Chapman-Enskog expansion method is 
found to be compatible with the underlying Boltzmann equation. 

Using Eq.'s (4-1) and (4-44), we can rewrite j££ = j[°^ a + 5J^ t as 

= e it" u v - {p - C V ■ u) A» v + 2 V A^ pa V p u a , (4-68) 
N< x = nu> x + X J-V^. (4-69) 

We should notice that T^ u and in Eq.'s (4-68) and (4-69) completely agree with 
the energy-frame currents proposed by Landau and Lifshitz. 12 ) Indeed, the respective 
dissipative parts 5T^ U and 5N^ in Eq.'s (4-68) and (4-69) meet Landau-Lifshitz's 
constraints given by Eq.'s (IT), (1-2), and (1-4). This was actually anticipated: In 
fact, if we take the simplest choice = u^, Eq.(2-31) tells us that the physical 
quantity transported by each particle in the system reads 

p ■ u. (4-70) 

To make clearer the physical meaning of (p ■ u) , we take the non-relativistic limit of 
this quantity: 



m 

p • u ~ m + — 



P 

u 

m 



(4-71) 



where = (u°, u) and p^ = p). Equation (4-71) shows that {p • u) can be 
identified as the kinetic energy of the fluid component measured in the rest frame of 
■u M . Thus, it is natural that the currents Jf s " with the choice ap 1 = becomes ones 
in the energy frame adopted by Landau and Lifshitz. 

4.3. Particle frame 

Another simple choice 9 = tt/2, i.e., ap 1 = m/(p ■ u)u^, gives the following 
constitutive equations 

be = -3C(37-4)~ 2 V-u, (4-72) 
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<5p = _£(3 7 -4)- 2 V -u, (4-73) 

Q^" = \ V^T, (4-74) 

77^ = 2r/Z\^ CT V p n (7 , (4-75) 

m<5n = 0, (4-76) 

mv» = 0. (4-77) 

We have utilized the following relations 

F 1 ( 7 r/2) = -3(3 7 -4)- 1 , (4-78) 

F 2 (7r/2) = -(3 7 -4)- 1 , (4-79) 

F 3 (tt/2) = 1, (4-80) 

Zi(tt/2) = 0, (4-81) 

Z 2 (tt/2) = 0, (4-82) 

which are derived from Eq.'s (4- 16)- (4- 20), respectively. We see that = while 
Q M / 0, which tells us that our generic equation with 9 = tt/2 becomes the hy- 
drodynamic equation in the particle frame or Eckart frame. We also note that the 
dissipative internal energy be and pressure 5p are finite in this frame in our formal- 
ism while Sn = 0; we remark that the presence of 5p was also the case in the energy 
frame. 

In terms of T^ v and iV 74 , we have 

T"" = (e-3C(37-4)" 2 V-u) vPu v - (p - ( (3 7 - 4)" 2 V • uj A» v 

+ A K VT + u u V^T) + 277 ^ vpu V p u a , (4-83) 
N^ = nu^. (4-84) 

It is noteworthy that the set of Eq.'s (4-83) and (4-84) is different from the phe- 
nomenological equations in the particle frame by Eckart, but may be regarded as a 
corrected version of that by Stewart, 14 ) which is derived from the relativistic Boltz- 
mann equation on the basis of the moment method: Indeed, the respective dissipative 
parts 5T^ U and <5iV> in Eq.'s (4-83) and (4-84) satisfy the conditions (1-5), (1-2), and 
(1-3) as Stewart's equation does, but not Eckart's constraints given by Eq.'s (1-1), 
(1-2), and (1-3). The latter also implies that the conditions of fit (2-39) imposed in 
the usual Chapman-Enskog expansion method is not compatible with the underlying 
Boltzmann equation. 

Now it is interesting that the above constraints on the energy-momentum tensor 
do not necessarily determine the detailed structure of it. In fact, Eq.(4-83) is different 
from that of Stewart, which is obtained from our equation through the replacements 

-C (37 - 4)" 2 V • u -> +C (37 - 4) _1 V • u, (4-85) 
A V^T — >■ A (V M T — T Du^), (4-86) 

with D = u^ 8^, although both of them satisfy the constraints (1-5), (1-2), and (1-3). 
There are two noteworthy points in Eq.(4-83): First, the sign of the thermodynamic 
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force with the bulk viscosity in our equation is the same as that in the Landau- 
Lifshitz equation and opposite that in the Stewart equation, because 7 > 4/3. ^ 
Second, the heat flow A^ u u p 5T vp = Q M does not contain the time-derivative term, 
Du p , i.e., the acceleration term, which is often interpreted as one of the important 
relativistic effects. u ) 

Now a natural question is which equation is preferable, ours or Stewart's or 
Eckart's? It is known that there have been no hydrodynamic equation including 
Eckart's and Stewart's which has the stable thermal equilibrium state, as is shown 
by Lindblum and Hiskock. 17 ) Owing to the absence of Du p , we can show that 
the thermal equilibrium solution of our equation becomes stable against a small 
perturbation, 37 ) in contrast to the Eckart and the original Stewart equations. Thus, 
we claim that our equation (4-83) is the proper hydrodynamic equation free from the 
instability problem in the particle frame. We emphasize that such an equation has 
been obtained for the first time, as far as we are aware of. A detailed and general 
account of this stability issue will be given in §5.3. 

4.4. A new frame 

The other interesting choice is 9 = 7r/4, for which we have ap 1 = {{jp ■ u) + 
m) I (y/2p ■ -u) and 

be = 0, (4-87) 

Sp = 0, (4-88) 

Q» = \(^^) 2 (Vr - z^TV^, (4-89) 
^h + z' V TJ 

= 2-qA» vpa V p u a , (4-90) 
mSn = 0, (4-91) 

m / = _A(A) 2 fv f 'T-z- 1 TV3. (4-92) 
+ z' ^ TJ 

Here, we have used the following relations obtained from Eq.'s (4- 16)- (4- 20), 



^(tt/4) = s -, (4-93) 

z 2 + z 2 (3 7 - 4) - 3 z [1 - (h - 1) (7 - 1)] 

F 2 (7r/4) = 0, (4-94) 

Y 3 („/4) = -J-, (4-95) 
h + z 

Zi(tt/4) = — r s -, (4-96) 

z 2 + z 2 (3 7 - 4) - 3z [1 - (h - 1) (7 - 1)] 

Z 2 (ir/4) = (4-97) 
h + z 

We note that the present choice 9 = 7r/4 gives a local rest frame where Q p 7^ 
and 7^ 0; recall that Q p = {y^ = 0) in the energy (particle) frame. Thus, our 
generic hydrodynamic equation with 9 = 7r/4 gives the equation in a local rest frame 
which is neither the energy frame nor the particle frame. We also remark that the 
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dissipative internal energy, pressure, and particle-number density are absent in this 
frame; 5e = 5p = 5n = 0. 

The alternative expressions in terms of T^ v and read 

T» v = eu^u v -pA» v 

+ 2r]A^ p<T V p u a , (4-98) 

ATM = nu » _ A ( ^—Y ( VT - z" 1 T V^) . (4-99) 
^h + z' V 

We see another unique feature of this frame that the bulk pressure term is absent, 
and hence £ does not play any role in this frame. As far as we know, this type of the 
relativistic dissipative hydrodynamics is written down for the first time, which was 
made possible by the introduction of the macroscopic-frame vector in the powerful 
RG method *). 

§5. Discussions 

In this section, we shall examine the properties of the derived equation in the 
generic frame and its advantageous nature over the existing equations in the liter- 
ature, focusing on transport coefficients and local rest frames. We give a general 
proof without recourse to numerical calculations that the hydrodynamic equations 
obtained in our formalism have the stable equilibrium state on the basis of the pos- 
itive definiteness of the inner product. 

5.1. Frame independence of transport coefficients 

In this subsection, we shall show that the transport coefficients (, A, and rj in 
Eq.'s (4-27)-(4-29) obtained in our formalism are independent of 9. For this purpose, 
we introduce a new linearized collision operator 

= ~ 2! 2 Zp S i Yl 3 > , Pa) C {S Pq + S piq - 5 P2q - S P3q ), 

Pi ^! p 2 2 Pa 3 

(5-1) 

which is independent of 9. We note that the inverse matrix of C pq reads 

££=L£(q-e q )-\ (5-2) 

Using £~ 1 and the ^-independent microscopic representations of the dissipative cur- 
rents (n p , Jp, TTp U ) in Eq.'s(4-13)-(4-15), we can represent the transport coefficients 

In a previous short communication, 36 ' we discussed the equation obtained with an erroneous 
setting 9 = — it /A. We note that this choice is not appropriate and the resulting equation should 
be abandoned since — 7r/4 is out of range of 8 (4-7) that guarantees the positive definiteness of the 
inner product. 
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as 

( = -^(n,c- 1 n) cq , (5-3) 

A = ^(J^,£- 1 J M ) cq , (5-4) 
V = -^ f (^,C- 1 n fll/ ) cq , (5-5) 
where we have introduced the new inner product 

v 

It is noted that this inner product is identically the thermal average, which is inde- 
pendent of 6. Equations (5-3)-(5-5) show that £, A, and r] are ^-invariant, and agree 
with the definitions used in the literature. 

5.2. Absence of macroscopic-frame vector leading to Eckart's particle-frame equa- 
tion 

In §4, we have seen that the hydrodynamic equation in the particle frame does 
not take the form of Eckart, if the equation is to be derived from the relativistic Boltz- 
mann equation. This fact might not be well known apart from some exceptions. 15 - ) 
For completeness and instructive purpose, we shall here give a proof that there can 
not exist a macroscopic-frame vector ap 1 that leads to Eckart's constraints, (IT), 
(1-2), and (1-3). It will be found in the course of the proof that the orthogonality 
condition between the Po and Qo space exactly corresponds to the phenomenologi- 
cal constraints imposed to the dissipative part of the energy-momentum tensor 8T^ V 
and the particle current 5N^ defined in Eq.'s (3-59) and (3-61). 

We first note that the dissipative currents SJ^ given in Eq.(3-54) can be ex- 
pressed as 

<st=E^^o P / P cq P , (5-7) 
p 

where 4> p is defined by 

^^-[L^Qof^F]^ [rVr],^- (5-8) 

where ip%£ are the first-excited modes (3-67) and X^ a the thermodynamic forces 
defined by Eq.(3-63). Note that 4> p belongs to the Qo space and is accordingly 
orthogonal to the five vectors ifQ p belonging to the Po space; 

(pg, 4>) =0 for a = 0, 1, 2, 3, 4. (5-9) 

Here, the inner product is defined by Eq.(3T7) where ap 1 enters in the form of (p-a p ). 
We are going to show that these five identical equations exactly correspond to the 
ansatz's on 5T^ U and SN^; see Eq.'s (3-59) and (3-61). 
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Firstly, let us take the case of the energy frame with the choice (A) a p = u^. 
Then, Eq.(5-9) is reduced to 

E \{p-u)r*^ P 4>v = ^ (s-io) 

V 

which actually expresses a set of equations as follows 

u u 5T» v = 0, (5-11) 
Up 5N^ = 0. (5-12) 

Note that Eq.(5Tl) implies the following two equations, 

Se = Ufl u u ST^ = 0, (5-13) 
Q p = A^u u ST^ = 0. (5-14) 

Thus, one can readily see that these equations coincide with Landau-Lifshitz's ansatz's, 
(1-1), (1-2), and (1-4). 

Similarly, with the choice of (B) a p = m/(p ■ u) u^, we have 



E 4)^/ P Cq ^ P ^ = 0, (5-15) 



p 

p 



which means that 



5N» = 0, (5-16) 
ST^ = 0. (5-17) 

Here, we have used the on-shell condition m 2 = p^ p p . Equation (5-16) is equivalent 
to the set of equations 

Sn = u p 5NV = 0, (5-18) 

u u = A^SN" = 0. (5-19) 

Thus, one sees that Stewart's ansatz's, (1-5), (1-2), and (1-3), are derived in the 
particle frame. 

It is now easy to see that there exists no macroscopic-frame vector a p leading 
to Eckart's ansatz's given by Eq.'s (1-1), (1-2), and (1-3), simultaneously. Indeed, in 
order to lead to the ansatz's (1-2) and (1-3) on 5N^, Eq.(5-9) with a = fi requires 
that (p • a p ) is independent of i.e., 

(p ■ a p ) = const.. (5-20) 

On the other hand, in order to lead to the ansatz (1-1) on 5T^ V , Eq.(5-9) with a = 4 
must lead to 

(p ■ a p ) = const, x (p ■ u) 2 , (5-21) 
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which is in contradiction with Eq.(5-20). Thus, we conclude that there exists no a,p 
leading to Eckart's ansatz's (l-l)-(l-3) simultaneously. 

In short, the phenomenological ansatz's on the dissipative parts of the energy- 
momentum tensor ST^ V and the particle current SN^ have a definite correspondence 
to the orthogonality conditions posed to the dissipative part of the distribution 
function as a solution to the underlying Boltzmann equation. In particular, whenever 
the particle frame is taken where the particle flow is constructed so as to have no 
dissipative part, the dissipative part of the energy- momentum tensor must satisfy 
Eq.(l-5) but not Eq.(l-l). This fact implies that the Eckart equation can not be the 
hydrodynamic equation in the particle frame as the slow (or infrared) dynamics of 
the relativistic Boltzmann equation, and hence nor have microscopic foundation. 

5.3. Stability analysis of steady solution of relativistic hydrodynamic equation 

In this subsection, we shall give a general proof that steady solutions of our 
generic relativistic dissipative hydrodynamic equation (3-68) are stable against a 
small perturbation, on account of the positive definiteness of the inner product (3-18). 
Here, a steady solution means that it describes a system having a finite homogeneous 
flow with a constant temperature and a constant chemical potential, as follows, 

T(a-T)=T , (5-22) 
h(<t;t) = hq, (5-23) 
U/x(o- ; r) = u ^, (5-24) 

where To, no, and uq^ are constant. We note that the steady states include the 
thermal equilibrium state as a special case. We remark that such a stability of 
our equation was demonstrated in a previous paper by the present authors 37 ) by 
a numerical calculation for a rarefied gas with some specific models for the cross 
section in the collision integral. The general proof of the stability of our equation is 
given for the first time in the present paper. 

To show the stability of the steady solution, we apply the so-called linear stability 
analysis to the relativistic dissipative hydrodynamic equation (3-68). We expand T, 
fi and around the steady solution as follows: 

T(a; t) = T + ST(a; r), (5-25) 
fi(a; t) = fi + 5fi(a; t), (5-26) 
Up(a; t) =u p + 5up(a; t). (5-27) 

We assume that the higher-order terms than the second order with respect to ST, 
Sfj,, and 5up can be neglected since these quantities are small. For convenience, we 
introduce the new variables SX a composed of ST, Sfi, and Su^ by 

( -S(ua/T) = -5up/T + 5Tuqp/Tq for a = fi, 

SX a = { / , „ 5-28 

I m- 1 S(fi/T) = m" 1 {Sn/T Q - ST // /T 2 ) for a = 4. 



Substituting Eq.(5-28) into Eq.(3-68), we obtain the linearized equation governing 
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5X a as 

+ (( iff* , 4 ) V M + ( (ft* , L- 1 <p? > V, V u ) SXp = 0, (5-29) 

with iflp = p fJ, (fQ p /(p ■ a p ) defined in Eq.(3-55). Here, we have used the following 
simple relations 

S(f?) = f?<P% p 6X a , (5-30) 
5{X m ) = V M <5X a . (5-31) 

We note that all of the coefficients in Eq.(5-29) take a value of the steady solution 
(T, fx, Uu) = (To, /Uo, uq^). With use of the orthogonality condition between the 
Po and Qo spaces and the definitions of rf^f and rff au ^ in Eq.'s (3-26) and (3-66), 
respectively, we can reduce Eq.(5-29) to 



A ' 13 -^$Xp + B a ' p 5X(s = 0, (5-32) 



where A a, @ and B a, @ are defined by 

A^^r,f, (5-33) 
B^ ee < ff a , 4 ) V„ + V M V u . (5-34) 

We convert Eq.(5-32) into the algebraic equation, using the Fourier and Laplace 
transformations with respect to the spatial variable and the temporal variable r, 
respectively: By inserting 

5X a (a ; t) = 6X a (k ■ A) e ik — A ^ (5.35) 

into Eq.(5-32), we have 

(AA a >P -B a >P)8Xp = 0, (5-36) 

where B a ^ is defined by 

B a * = i ( ff a , 4 ) k, - r,r P K K- (5-37) 

In the rest of this section, we use the matrix representation when no misunderstand- 
ing is expected. Since we are interested in a nonvanishing solution SX / 0, we 
impose 

det(AA- B) = 0, (5-38) 
which leads to the dispersion relation 

A = A(k). (5-39) 
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The stability of the steady solution (5-22)-(5-24) against a small perturbation means 
that SX will not increase under time evolution. Therefore, our next task is to show 
that the real part of A(k) is non-negative for any k^ 1 . 

Now we first note that the matrix A is a real symmetric and positive-definite 
matrix: 

= ( if , ip ) > for w a ^ 0, (5-40) 

with (p p = w a <Pq . Here, we have used the positive definiteness of the inner product 
(3T8). Equation (5-40) means that the inverse matrix A~ l exists, and A~ l is also 
a real symmetric positive-definite matrix. Thus, using the Cholesky decomposition, 
we can represent A' 1 as 

A~ 1 = t UU, (5-41) 

where U denotes a real matrix and t U a transposed matrix of U. Substituting 
Eq.(5-41) into Eq.(5-38), we have 

det{A-UB t U) = 0. (5-42) 

It is noted that A(k) is an eigen value of U B l U. 

We notice the following theorem: The real part of the eigen value of a complex 
matrix C is non-negative when the hermite matrix Re(C) = (C + C')/2 is semi- 
positive definite. Applying this theorem to the present case, we find that the real 
part of A(k) becomes non-negative for any k^ 1 when Re(U B t U) is a semi-positive 
definite matrix. In fact, we can show that Re(U -B*C7) is semi-positive definite as 
follows: 

w a [Re(U B t U)] a ^ wp = w a [U Re(B) l U] a ^ wp 

= [w U] a [Re(B)] a,/3 [wU] p 
= -[w U] a rj{ avfi k^ k v [w U]p 
= -( kfj, [w U] a ip1 a , L^ 1 k v [w U]p iff ) 
= -(V, L" 1 rj)) > for w a / 0, (5-43) 

with ip p = k^ [w U] a ip^p . Therefore, we conclude that the steady solution in Eq.'s 
(5-22)-(5-24) is stable against a small perturbation. 

We now see that the relativistic dissipative hydrodynamic equation (3-68) ob- 
tained by the RG method has a stable steady solution in a definite way. It is note- 
worthy that this property is kept in the several equations derived in the setting of 
ap 1 = 9p, i.e., the energy-frame equation in Eq.'s (4-68) and (4-69), the particle-frame 
equations in Eq.'s (4-83) and (4-84), and so on. We stress that this is for the first 
time that the relativistic dissipative hydrodynamic equation in the particle frame 
has been obtained whose steady solution is stable *) . 

*-* In the previous rapid communication, 37 ' we discussed the stability only of the thermal- 
equilibrium solution, i.e., u 0fl = (1, 0, 0, 0), with use of a rarefied-gas approximation, where a 
differential cross section in the collision integral is treated as a constant. 
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§6. Summary and Concluding Remarks 

In this paper, we have given a full and detailed account of the derivation of the 
first-order relativistic dissipative hydrodynamic equations as a reduction of dynamics 
from the relativistic Boltzmann equation with no heuristic assumptions, such as the 
so-called conditions of fit used in the standard methods. This was made possible by 
adopting a powerful reduction theory of dynamics, i.e., the renormalization-group 
method and by introducing the macroscopic-frame vector which defines the macro- 
scopic local rest frame; thereby we successfully have a coarse-grained and Lorentz- 
covariant equations in a generic frame. 

The five hydrodynamic modes are naturally identified with the same number of 
the zero modes of the linearized collision operator. The excited modes which are 
to modify the local equilibrium distribution function and give the dissipative terms 
are defined so that they are precisely orthogonal to the zero modes with a properly 
defined inner product for the distribution functions. It is worth emphasizing that 
the dissipative terms are constructed without recourse to the so-called conditions 
of fit which are imposed to the dissipative terms in an ad-hoc way in the standard 
methods. In our method, the validity of the conditions of fit is checked as a property 
of the derived equations. On the basis of the nice properties of the properly defined 
inner product and the very nature of the hydrodynamic modes as the zero modes 
of the linearized collision operator, we have shown that the so-called Burnett term 
does not affect the hydrodynamic equations. 

Our hydrodynamic equation with the dissipative currents given by Eq.(4-44), 
which still contains the macroscopic-frame vector, is a master equation for a generic 
local rest frame, and derives hydrodynamic equations in various local rest frames 
with a specific choice of the macroscopic- frame vector. We have shown that our 
energy- momentum tensor and particle current in the energy frame, given by Eq.'s 
(4-68) and (4-69), respectively, coincide with those proposed by Landau and Lifshitz, 
but those in the particle frame, given by Eq.'s (4-83) and (4-84), are slightly different 
from those given by Eckart and by Stewart. Our generic equation is also reduced to 
a novel relativistic hydrodynamic equation with the energy-momentum tensor and 
the particle current given by Eq.'s (4-98) and (4-98), respectively, where the bulk 
pressure (viscosity) term is absent. 

Furthermore, we have proved that the derived equation in a generic local rest 
frame has a stable equilibrium state owing to the positive definiteness of the inner 
product, although such a generic stability was suggested by a numerical calculation 
for some specific parameter sets in a previous paper by the present authors. It is 
worth emphasizing that all of our equations have a stable equilibrium state even in 
the particle frame in contrast to the Eckart and Stewart equations; note that such 
a drawback of the equation that the thermal equilibrium state can be unstable is 
taken over to some causal equations including the Israel-Stewart equation. 

In conformity with the above fact, we have proved that the Eckart equation can 
not be compatible with the underlying relativistic Boltzmann equation. This proof 
is based on the following significant observation that the orthogonality condition 
of the excited modes to the zero modes coincides with the ansatz's posed on the 



38 



K. Tsumura and T. Kunihiro 



dissipative parts of the energy-momentum tensor and the particle current in the 
phenomenological equations. 

Thus, we have arrived at a sound starting point for attacking the last problem, 
i.e., constructing a causal relativistic hydrodynamic equation: We can expect that 
the present method applies to derive a causal equation on a sound basis, without re- 
course to any ansatz's as is done in other method such as the Maxwell-Grad moment 
method. 18 )' 20 ^' 26 ) We have found 21 ) that a non-trivial extension of the Po space in 
our formalism gives a causal equation, whose low-frequency limit is identical to the 
equation derived in this paper. This extension to derive causal equations will be 
reported in the forthcoming paper, 21 ) where a correct moment method will be also 
proposed. 

We stress that all of the equations derived in this paper are consistent with the 
underlying kinetic equation, so the equations and also the method developed here 
may be useful for the analysis of the system where the proper dynamics describing 
the system changes from the hydrodynamic to kinetic regime, as in the system near 
the freeze-out region in the RHIC phenomenology. 42 )~ 44 ) 

Due to the absence of the problematic Burnett term in our relativistic equations, 
we expect that applying the RG method to the non-relativistic Boltzmann equation 
leads to the Navier-Stokes equation without the Burnett term, whose existence is 
known to be inevitable 39 ) when the Chapman-Enskog expansion method is applied. 

Finally, we emphasize that our method itself has a universal nature and can 
be applied to derive a slow dynamics from kinetic equations other than the simple 
Boltzmann equation. For example, it is successfully applicable to derive the dissi- 
pative hydrodynamic equation for the multi-component system in the energy frame 
from the multi-component Boltzmann equation. 18 ) Furthermore, it would be inter- 
esting to apply the present method for extracting a Lorentz-covariant hydrodynam- 
ics of strongly interacting systems out of the equilibrium from the Kadanoff-Baym 
quantum-transport equation 45 ) for the non-equilibrium many-body Green's function. 
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Appendix A 

Most general form of phenomenological relativistic dissipative hydrodynamic 
equations with dissipative internal energy 5e, pressure dp, and density 5n 



The most general form of the dissipative parts of the energy-momentum tensor 
and particle current takes the following forms, 

8T» U = 5e U v U u -5pA» v + Q»u u + Q v u» + nv u , (A-l) 
SN» = Snu» + ^, (A-2) 

respectively. Note that we have made it explicit that the dissipations may cause 
corrections to the internal energy Se, the pressure Sp and the particle-number den- 
sity Sn, although only 5p is usually considered (as the bulk pressure term) in the 
literature. We remark that Se and Sp are related to the dissipative part of the 
energy-momentum tensor, as follows, 

u a ST» V u v = Se, (A-3) 
ST^ = 5e-3Sp. (A-4) 

In the present parametrization, e + Se, p + Sp, and n + Sn are the internal 
energy, pressure, and particle-number density in the dissipative system, with e = 
e(T, fi), p = p(T, fi), and n = n(T, fi) being the corresponding quantities in the 
local equilibrium state with the temperature T and the chemical potential fi. In the 
tensor decomposition, 

Q» = A^uPT up , (A-5) 
V » = A^N V , (A-6) 
n^ = A^ a T pa , (A-7) 

where A^ upa is the space- like, symmetric and traceless tensor defined in Eq.(4-23). 

Owing to the properties of the dissipative parts, Q^u a = 0, u a = 0, II ^ = 
n vp , and II v = U^p, = 0, the total number of independent components of Q^, u^, 
and n pv is eleven. Since T pv and have fourteen components in total, Se, Sp, and 
Sn can not be independent of each other, but the number of independent component 
is one other than T and /i. We take Sp = II as the independent component as a 
natural choice, then Se and Sn may be expressed as Se = f e H and Sn = f n II, where 
f e and f n are functions of T and [i; f e = f e (T, fi) and /„ = / n (T, /i). Here, we have 
assumed that the dissipative order of Se and Sn is the same as that of Sp at most. 
We remark that in terms of f e , 

u a ST^u u = f e II, (A-8) 
ST^ = (f e -3)n. (A9) 

Now we shall show that the usual phenomenological derivation of the hydrody- 
namic equations allows the existence of Se and Sn, i.e., finite values of f e and /„, in 
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the relativistic dissipative hydro-dynamic equations. We emphasize that it can not 
be excluded on the general ground that f e and f n may have finite values, although 
the phenomenological theory can not determine their parametric forms as functions 
of T and ji. Although the first-order equation is considered in this article, the same 
argument equally applies to the second-order equations. 
Now the entropy current is given by 



(A-10) 



The second law of thermodynamics reads d^S^ > 0. 
The divergence of is found to take the form 



= n 



f e D±-±V^-f n D^ 



- D Ufl + 



(A-ll) 



where D = u^ and V M = A^ v d u . Here, we have used the hydrodynamic equation 
and the first law of thermodynamics, 



D(p/T) + eD(l/T) -nD(ji/T) = 0. 
In the particle frame where = 0, Eq.(A-ll) is reduced to 



(A-12) 



dpS" = n 



f e D±-±V» UfM -f n D^ 



1 



+ 



(A-13) 



For assuring the second law of thermodynamics, one should make the right-hand 
side semi-positive definite. A natural choice is the following constitutive equations, 



77 = CT 



f e D±-±Wu„-f n D^ 



= -XT 2 



T T 



(A-14) 

(A-15) 
(A-16) 



with £, A, and r] being the bulk viscosity, heat conductivity, and shear viscosity, 
respectively. Indeed, the above constitutive equations lead to 



CT XT 2 2 V T 



> 0. 



(A-17) 



Thus, one sees that the relativistic dissipative hydrodynamic equations with 
finite f e and f n , or equivalently with finite Se and Sn, is compatible with the second 
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law of thermodynamics. With a special choice, f e = fn = 0, Eq.'s (AT4)-(AT6) lead 
to the constitutive equations proposed by Eckart that are commonly used. 

In the case of the energy frame where = 0, a similar argument leads to Eq.'s 
(A-14), (A-16), and 

= A/r 2 V^, (A- 18) 

with h = {e + p)/nT being the reduced enthalpy per particle. We remark that these 
equations are reduced to the constitutive equations by Landau and Lifshitz only 
when one can set f e = f n = 0. 

The phenomenological theory can not proceed further to specify f e and f n ; 
the values of f e and f n can be determined only from a microscopic theory. In the 
present work, we show that the microscopic theory leads to f e = f n = in the energy 
frame, while f e = 3 and f n = in the particle frame, implying that ST^ = but 
u^ST^Uv = 377^0. 



Appendix B 

Brief account of the renormalization- group method with an example 

In this Appendix, we first show how the renormalization-group (RG) method 28 )- 33 ) 
works as a reduction theory of the dynamics, adopting van der Pol equation with a 
limit cycle, as an example. Then, we briefly present a foundation of the RG method. 
A detailed account of the RG method including its foundation may be seen in Ref.'s 
29), 30), 32), 33) or review articles 35), 46). 

B.l. RG method applied to van der Pol equation 

Let us take the van der Pol equation which admits a limit cycle: 

x + x = e(l - x 2 )x, (B-l) 

where e is supposed to be small. 

Let us solve the same problem by the RG method. Let x(t; to) De a local solution 
around t ~ Vto, and represent it as a perturbation series; x(t; to) = Xo(t; to) + 
exi(t; to) + e 2 X2(t; to) + ■ ■ ■ . In the RG method, the initial value W(to) matters: 
We suppose that an exact solution is given by x(t) and the initial value of x(t; to) at 
t = to is set up to be x(to); i.e., W(to) = x(to; to) = x(to)- The initial value as the 

exact solution should be also expanded as W(to) = Wo(t )+e H^i(to)+e 2 W^f^H • 

The RG method is actually the method to obtain the initial value W(t) as an exact 
solution or approximate solution valid in a global domain asymptotically. 

The zeroth-order equation reads 



Cxo = 

the solution to which can be written as 



2 



xo = 0, (B-2) 



x (t; t ) = A(t ) cos(t + 0(to)). 



(B-3) 
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Here, we have made it explicit that the integral constants A and 9 may depend on 
the initial time to- The equation for x\ reads 

( A 2 \ A 3 
Cxi = -A ( 1 - — J sin + — sin3<£(t), (B-4) 

with <fi(t) = t + 9o(to). Notice that the first term in the right-hand side is a zero 
mode of the linear operator £ appearing in the left-hand side. Thus, the special 
solution to this equation necessarily contains a secular term which is given by f 
times a zero mode of C. Since we have supposed that the initial value at t = to is on 
an exact solution, the corrections from the zeroth-order solution should be as small 
as possible. This condition is realized by setting the secular terms appearing in the 
higher orders vanish at t = to, which is possible because we can add freely zero mode 
solutions to a special solution. Thus, the first-order solution is uniquely written as 

xi(t;t ) = (t - t ) ^ (l-^j S m <f>{t) - ^ sin3 0(t). (B-5) 

Notice that the secular term surely vanishes at t = to in Eq.(B-5), implying that its 
initial value at t = to reads 

xi(t ;to) = -^^sin3^ ). (B-6) 

The perturbative solution up to this order reads x = xo + ex±, which becomes, 
however, invalid when \t — to I ~~ > large, because of the secular term. 

Now notice that the function x(t; to) corresponds to a curve drawn in the (t, x) 
plane for each to; in other words, we have a family of curves represented by x(t; to) in 
the (t, x) plane, a member of which is parametrized by to- An important observation 
is that each curve is close to the exact solution in the neighborhood of t = to- Thus, 
an idea is that the envelope curve of the family of curves should give a global solution. 
The envelope curve can be constructed by solving the following equation 



dx 

dto~ 



= 0, (B-7) 

to=t 



which leads to the following equations for A(t) and 4>(t), 



A = ei(l-^-), (B-8) 



2 V 4 J 

<p = l. (B-9) 

These equation are readily solved, and one sees that as t — > oo, A(t) — > 2 asymptot- 
ically, meaning the existence of a limit cycle with a radius 2. 

The resultant envelope function as a global solution is given by 

x E (t) = x(t; t) = W(t) = A(t) cos(t + O ) - e sin(3 1 + 3 0„), (B-10) 
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with A{t) being the solution of Eq.(B-8). Thus, we have succeeded in not only 
obtaining the asymptotic solution as whole but also extracting the slow variables 
A(t) and (p(t) explicitly and their governing equations. 

One now sees that when there exist zero modes of the unperturbed operator, 
the higher-order corrections may cause secular terms, which are renormalized into 
the integral constants in the zeroth-order solution by the RG/envelope equation 
(B-7), and thereby the would-be integral constants are lifted to dynamical but slow 
variables. 

In the derivation of hydrodynamic equations from the kinetic equation, the 
would-be integral constants corresponding to A and 4> are the temperature T, the 
chemical potential fi, and the flow velocity {vP = 1) characterizing the local 
equilibrium state, which are to be lifted to the slow dynamical variables and their 
governing equations are identically the hydrodynamic equation. 

B.2. Foundation of RG method 

Now we present a foundation to the RG method using a Wilsonian equation 47 ) 
or flow equation by Wegener. 48 ) 

Let us take the following n-dimensional equation; 

^=F(X,t), (B-ll) 

where n may be infinity. Let X(t) = W(t) be an yet unknown exact solution to 
Eq.(BTl), and we try to solve the equation with the initial condition at t = Vto; 

X(t = t ) = W(t ). (B-12) 

Then, the solution may be written as X(t; to, W(to)). 

Now the basis of the RG method lies in the fact that W(to) can be determined 
on the basis of a simple fact of differential equations. We notice that when the initial 
point is shifted to t' , the resultant solution should be the same as long as W(t) is 
an exact solution, i.e., 

X(t; t , W(t )) = X(t; t' , W(t' )). (B-13) 
Taking the limit t' Q — > to, we have 

dX dX OX dW „ 

= 1 = 0. (B-14) 

dt dt dW dto K ' 

This equation gives an evolution equation or the flow equation of the initial value 
W(to). This equation has the same form as and corresponds to the non-perturbative 
RG equations (flow equations) by Wilson, 47 ) Wegner-Houghton, 48 ) and so on, in 
quantum field theory and statistical physics: The 'initial time' to corresponds to (the 
logarithm of) the renormalization point. We emphasize that the equation (BT4) is 
exact; we have no recourse to any perturbation theory so far. 

The problem is to construct the seed of the RG equation (B-14), i.e., X(t; to, W(to)). 
For that, let us take the perturbation theory. In this case, X(t;to,W(to)) and 
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X(t;t 0l W(t )) may be valid only for t ~ to an d t ~ t , which condition is nat- 
urally satisfied when to < t < t (or t < t < to) because the limit t' — > to is 
taken eventually. Thus, when a perturbative expansion is employed for constructing 
X(t;to, W(to)), it is necessary to demand 



dX 
dto~ 



dX 



t =t 



dtn 



+ 



dX dW 



t =t 



dW dtr 



= 0, 



(B-15) 



t =t 



which automatically implies the condition that to = t- 

We have already suggested an interpretation that this equation can be identified 
as a condition to construct an envelope of the curves represented by the unperturbed 
solutions with different initial times to's: 29 ^ When to is varied, X(t;to, W(to)) gives 
a family of curves with to being a parameter characterizing curves. Then, Eq.(BT5) 
is a condition to construct the envelope of the family of curves which are valid only 
locally around t ~ to- The envelope is given by X(t; to = t) = W(t), i.e, the initial 
value. 

Now we shall show that X(t; to = t) = W(t) satisfies the original equation (B-ll) 
in a global domain up to the order with which X(t; to) satisfies around t ~ to- Let 
X(t;to) is an approximate solution to Eq.(B-ll) around t ~ t ; 

dX(t;t ) 



Then, we have 



dt 



dW(t) dX(t;t c 



F(X(t;t ),t). 



(B-16) 



dt 



dt 

dX(t;t ) 



+ 



dX(t-t ) 



t =t 



dtn 



t =t 



dt 



to=t 

~F(X(t;to), t)\ t0=t , 
= F(W(t),t), 

on account of Eq.(B-15). This proves the above statement. 



(B-17) 



Appendix C 

Derivation of first- excited modes with a generic macroscopic-frame Vector 



In this Appendix, we explicitly derive the first-excited modes (p^ in the generic 
local rest frame with the macroscopic-frame vector being specified by Eq.(4-6), i.e., 
dp = Op = [((p ■ u) cosO + m sin Q)/(p ■ u)] u^. 

For convenience, we introduce a dimensionless quantity 18 ) dependent on z = 
m/T, 

v 
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This quantity for I = 0, 1, 2, 3, 4, 5 reads 

a = z~ 2 (a-4), 
a x = 1, 

CL2 = h — 1, 

03 = 3 a + z 2 , 



a 4 = (15 + z 2 )n + 2z 2 , 

a 5 = 6 (15 + z 2 ) h + z 2 (15 + z 2 ) 



(C-2) 
(C-3) 
(C-4) 
(C-5) 
(C-6) 
(C-7) 



respectively. Here, a is defined in Eq.(4-21) and .fQ(z) denotes the modified Bessel 
function, 



K t (z) 



2 e £\ 



dr (r 2 - z 2 ) 



2^-1/2 e -r 



(C-8) 



which satisfies the recurrence relation 



K e+1 (z) = K t _ 1 (z) + — K t (z). 



(C-9) 



First, the setting of = 6^ leads us to the metric matrix t/q given by 



t]q = n 1 



(03 cos 9 + z a 2 sin 6*) a M ti^ 



z/4 



4/i 

Vo 
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+ ((z 2 ai - a 3 ) cos 6» + z (z 2 a - a 2 ) sin 6») - Z\ MJy 

nT 2 z(d2 cos 6* + z a\ sin #) u' 1 , 
n T 2 z 2 (ai cos 6* + z ao sin 6*) . 



By a tedious but straightforward manipulation, we have 

(nT 2 )- 1 (a(9)u^u v + B{9) A^, 

(nT 2 y 1 C(9)u l „ 
(nT 2 )- 1 D{9), 



where 



A(0) = 



a\ cos 9 + ao z sin 9 



(C-10) 

(C-ll) 
(C-12) 



(C-13) 
(C-14) 
(C-15) 



(a 3 ai — a 2 ) cos 2 6* + (03 ao — a 2 ai) z sin cos + (a 2 ao — a 2 ) z 2 sin 2 6* ' 

(C-16) 

m = TjT, ~mr7i r~7i9~~ -t— m> (C-17) 



(z 2 ai — 03) cos # + (z 2 ao — a 2 ) z sin ' 



c(<0 = -- 



a 2 cos 6* + ai z sin 



z (03 ai — a 2 ,) cos 2 9 + (03 ao — a 2 ai) z sin 9 cos + (a 2 ao — a 2 ) z 2 sin 2 9 ' 
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(C-18) 



D{0) = 



03 cos 9 + a 2 z sin # 



z 2 (03 a\ — a 2 ) cos 2 9 + (03 ao — a 2 ai) 2 sin# cos# + (02 ao — a\) z 2 sin 2 

(C-19) 



Next, the inner product ( , ip^ ) is evaluated to be 

( <pl , (ft* ) = n T 2 (a 3 u a u» u b + {z 2 ai - 03) ^ (u° + u" A ba + u b A a ^ , 

(C-20) 

( <fQ , (p^ A ) = nT 2 (z a 2 u a u» + z (z 2 a - o 2 ) ^ , (C-21) 
( ^ , ) = n T 2 (z a 2 u»u b + z (z 2 a - a 2 ) i Z^ 6 ) , 



(C-22) 
(C-23) 



Finally, combining 77^ in Eq.'s (C-13)-(C-15) and (tp" , pf) in Eq.'s (C-20)- 
(C-23) together with 

= JT P ( p " ^ " {p " 6p) ^ * ( ^ ' ^ } ) ' (C ' 24) 



we have 



ua 



P-O* 



p-9 p 



n p (u»u' / Y 1 (9)-A^Y 2 (9^ 

+ (J£ u" + j; u") y 3 (0) + tt^ for a = n, (C-25) 
n p u^ Z x {9) + Z 2 {9) (01 n - i 



Here, we have introduced the following quantities 



(02 ao — a 2 ) z 2 (p • «) 2 — (03 ao — a 2 a\) zm(p • u) + (03 a\ — a 2 ) m 2 

-3 (a 3 ai - a 2 ,) 

A^ p u [z 2 ao — 02) z (p- u) — (z 2 a\ — a 3 ) m 



n„ = 



jp 



Yi(0) 



-{z 2 a - a 2 ) z 

-3 (a 3 ai - a 2 .) sin 2 6* 



(C-26) 

(C-27) 
(C-28) 



(03 ai — a 2 ) cos 2 + (03 ao — a 2 a±) z sin 9 cos + (a 2 ao — a 2 ) z 2 sin 2 ' 

(C-29) 
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Y 2 (9) = 



(03 a\ — a 2 ) (cos 2 9 — sin 2 



(03 ai — a 2 ) cos 2 9 + (03 ao — a 2 a\) z sin 9 cos # + (02 ao — a 2 ) z 2 sin 2 9 ' 

(C-30) 

v rfl N -(z 2 a -a 2 )z sinfl 

(z z a\ — 03) cos 9 + (z z ao — a 2 ) z sin & 



3 (03 ai — a 2 .) cos 9 sin 



(03 ai — a 2 ) cos 2 6* + (03 ao — a 2 a±) z sin 9 cos 9 + (a 2 ao — a 2 ) z 2 sin 2 6* ' 

(C-32) 

z m s r-^s - (oss) 

(z z ai — 03) cos + (z 2 ao — 02) 2 sm 9 

Substituting ao, a\, a 2 , and 03 in Eq.'s (C-2)-(C-7) into the above equations, we 
arrive at the explicit representation of 99^ shown in Eq.(4-ll). 
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